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ABSTRACT 



Stern et all ( 120121 ) presented a study of WISE selection of AGN in the 2 deg 2 



COSMOS field, finding that a simple criterion Wl-W2>0.8 provides a highly 
reliable and complete AGN sample for W2<15.05, where the Wl and W2 pass- 
bands are centered at 3.4/zm and 4.6/zm, respectively. Here we extend this study 
using the larger 9 deg 2 NOAO Deep Wide-Field Survey Bootes field which also 
has considerably deeper WISE observations than the COSMOS field, and find 
that this simple color-cut significantly loses reliability at fainter fluxes. We define 
a modified selection criterion combining the Wl— W2 color and the W2 magni- 
tude to provide highly reliable or highly complete AGN samples for fainter WISE 



1 Jet Propulsion Laboratory, California Institute of Technology, 4800 Oak Grove Drive, Mail Stop 169-236, 
Pasadena, CA 91109, USA. E-mail: roberto. j .assef@jpl.nasa.gov 

2 NASA Postdoctoral Program Fellow 

department of Astronomy, The Ohio State University, 140 W. 18th Ave., Columbus, OH 43210, USA 

4 The Center for Cosmology and Astroparticle Physics, The Ohio State University, 191 West Woodruff 
Avenue, Columbus, OH 43210, USA 

5 Physics & Astronomy, University of Leicester, 1 University Road, Leicester, LEI 7RH, UK 

6 Department of Physics, University of Missouri, 5110 Rockhill Road, Kansas City, MO 64110, USA 

7 School of Physics, Monash University, Clayton 3800, Victoria, Australia 

8 Infrared Processing and Analysis Center, California Institute of Technology, Pasadena, CA 91125, USA 
9 Steward Observatory, University of Arizona, Tucson, AZ 85721, USA 

10 Department of Physics, University of California, One Shields Avenue, Davis, CA 95616, USA 

n Institute of Geophysics and Planetary Physics, Lawrence Livermore National Laboratory, Livermore, 
CA 94550, USA 



-2- 



sources. In particular, we define a color-magnitude cut that finds 130±4 deg -2 
AGN candidates for W2<17.11 with 90% reliability. Using the extensive UV 
through mid-IR broad-band photometry available in this field, we study the spec- 
tral energy distributions of WISE AGN candidates. As expected, the WISE AGN 
selection is biased towards objects where the AGN dominates the bolometric lu- 
minosity output, and that it can identify highly obscured AGN. We study the 
distribution of reddening in the AGN sample and discuss a formalism to account 
fo r sample incompl e tenes s based on the step-wise maximum-likelihood method 
of Efstathiou et al. Jl988[ ). The resulting dust obscuration distribution s depend 



strong ly on AGN luminosity, consistent with the trend expected for a ISimpson 



f l2005f ) receding torus. At L AGN ~ 3 x 10 44 erg s" 1 , 29±7% of AGN are observed 
as Type 1, while at ~ 4 x 10 45 erg s _1 the fraction is 64±13%. The distribution 
of obscuration values suggests that dust in the torus is present as both a diffuse 
medium and in optically thick clouds. 

Subject headings: galaxies: active — methods: statistical — quasars: general 



Introduction 



Active galactic nuclei (AGN) have been proposed to play an important role in sev- 
eral aspects of galaxy evolution, such as quenching star-formation in their host galaxies by 
heat ing and/or mechanic ally pushing their gas reservoirs into the intergalactic medium (see, 



Hop kins et al.l 120051 ). preventing cooling flows at the center of galaxy clusters (see, e.g., 



Croton et al.ll2006l ). an d possibly by contributing significantly to the reionization of the Uni- 
verse at high redshift (IGlikman et al.ll2010l 120111 ). Efficiently identifying AGN in all states 
of accretion and obscuration and accurately understanding their properties and structure is 
a key step to understand how galaxies evolve with cosmic time. 

AGN are among the most luminous objects in the Universe. While the bulk of the 
soft X-ray through NIR luminosity is thought to be thermal ly produced near the cent ral 
super- massive black hole (SMBH) in an accretion disk (e.g., IShakura fc Sunyaevl 119731 ). a 
considerable fraction of this luminosity is absorbed by dust surrounding the SMBH on parsec 
scales, and re-radiated in the m id-IR. This dust is generally believed to be responsible for the 



Type 1/2 dichotomy (see, e.g.. lAntonucci 



199 3J). The dust distrib ution is typically thought 



to have a quasi-toroidal shape (see, e.g., lUrry fc Padovanil Il995l ). with a s cale height an d 
opening angle that depends on the luminosity of the accretion disk (see, e.g.. ISimpsonl 120051 ). 
Because the dust is hot, at temperatures reaching its sublimation limit (~ 1500K), the 
emission is enhanced and much redder than the Rayleigh- Jeans tail in the 3 to lOyum range 
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(see, e.g., IStern et al.ll2005l ). Hence, even if dust in the torus is blocking our line of sight 
towards the SMBH and absorbing the soft X-rays through NIR emission of the accretion 
disk, the distinctive mid-IR warm dust emission from the torus is still observable. 



The Wide-field Infrared Survey Explorer (WISE; Wright et al.ll2010h is a NASA satellite 
with a 40cm aperture that imaged the whole sky in four mid-IR bands, centered at 3.4, 4.6, 
12 and 22/im. We refer to these bands as Wl, W2, W3 and W4, respectively. The fully 
cryogenic WISE science mission started in January 2010 and ended in August of the same 
year; all of this data has been publicly available since March 2012. While not one of its main 
goals, WISE is well suited to studying AGN as its bands are sensitive to their characteristic 
warm dust emission and are little affected by the obscuration expected from eithe r the dust in 



the to r us or in the interstellar medium of the host galaxy. In a companion work. IStern et al. 



(120121 . iPaper J ) investigated the power of WISE to identify AGN based solely on the Wl- 
W2 mid-IR color by comparing to known AGN in the COSMOS field. The selection also 
necessarily entails a flux cut, which is relatively shallow given the low ecliptic latitude of 
the COSMOS field and thus lower WISE coverage dept h. Using a set of AGN in the field 
selected from the Spitzer I nfrared Array Cam era (IRAC; iFazio et al.l 120041 ) colors according 
to the criteria developed by lStern et al.l (120051 ) as a control sample, we determined that at the 
depth of the COSMOS field a very simple selection criterion of Wl— W2>0.8 and W2<15.05 
produced an AGN sample with a contamination of only 5% and recovered nearly 80% of the 
IRAC-selected AGN in the field to that WISE depth. The 61.9 ± 5.4 deg" 2 space density 
of these AGN is about three times higher than that of similar bolometric lu minosity Type 1 



AGN found at optical wavelengths by the Sloan Digital Sky Survey (SP SS: iRichards et al. 



20021 ) . The color of this criterion is similar to that proposed earlier by lAssef et al. 



jioioL 



Wl-W2>0.85) using mock WIS E data constructed from SED models calibrated in this 
wavelength range by Spitzer data. lAssef et al.l (120101 ) have shown that WISE and SDSS are 
sensitive to AGN of the same bolometric luminosities for z < 4, implying that the i ncreased 
census is due to the sensitivity of WISE to those objects obscured by dust. In fact, IPaper I 
showed that the distribution of the X-ray hardness ratios of the WISE-selected AGN are, as 
expected, consistent with a considerable number of dust-obscured AGN. 

Because of the orbit and observing strategy of WISE, the depth of a field depends 
strongly on ecliptic latitude. The 2 deg 2 COSMOS field, close to the ecliptic, is represen tative 
of the shallowest WISE fields. In this work we extend the wor k presented in IPaper t to the 
much larger, 9 deg 2 NOAO Deep Wide-Field Survey (NDWFS: ljannuzi fc Devill999h Bootes 
field, which has considerably deeper WISE observations due to its higher ecliptic latitude. 
In §5] we describe the data sets and spectral ener gy distri bution (SED) fitting models which 
we use to extend the AGN selection criteria of IPaper il to fainter WISE fluxes in £j3j In 
§H] we study the broad-band SEDs of the WISE AGN candidates, and assess the accuracy 
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with which we can estimate photometric redshifts for them. Finally, in £j5l we study the 
distribution of the obscuring dust in AGN and present a method to correct for sample 
incompleteness due to reddening. Throughout this work we assume a flat ACDM cosmology 
with H = 73 km s -1 , Qm — 0.3 and Q\ = 0.7. We refer to all magnitudes in their native 
photometric system, i.e., AB for ugriz, FUV and NUV, and AO (Vega) for all other bands. 



2. Data and Modeling 



2.1. The NOAO Deep Wide-Field Survey Bootes Field 

NDWFS is a deep imaging survey in Bw, R, I and K of two 9 deg 2 fields in the 
constellations of Cetus and Bootes. We focus here on the Bootes field, for which follow-up 
deep imaging has been obtained for a wide range of wavelengths. Bootes also has deep and 
extensive spectroscopy. 

Follow-up ima ging of the Bootes field exists from the X-rays with Chandra (XBootes; 



Murray et al.ll2005f) to the ra dio from the Faint Images of the Radio Sky at Twe nty-centimeters 



(FIRST; iBecker et al.l |l995|) survey, the NRAO VLA S ky Survey (NVSS: ICondon et al. 
1998 ), the Westerbo rk Northern Sky Survey (WENSS; iRengelink et al.l Il997| ) and from 



de Vries et al.l ( 120021 ). The whole fi eld was observed with 90 seconds of exposure per po- 
sition in the IRA C Shallow Survey ( jEisenhardt et al. 2004 ). The Spitzer Deep, Wide-Field 
Survey fSDWFS: lAshby et al.l 120091 ) quadrupled this exposure, reaching 5a depths of 19.3, 
18.5, 16.3 and 15.6 mag for [3.6], [4.5], [5.8 ] and [8.0], respectiv ely. Additionally, we also use 
the Galaxy Evolution Explorer (GALEX; iMartin et al.l 120051 ) Deep Imaging Survey (DIS) 
a nd All - sky Im aging; Survey (AIS) FUV and NUV observations of the field, the z-band data 
of lCooll J2007h . the near-IR J, H and K s observations of NEWFIRM jGonzalez et al-lkoioh 
and the MIPS 24mm observations of the MIPS AGN and Galaxy Evolution Survey (MAGES; 
Jannuzi et al.ll2010l ). For our work, we use 6" aperture magnitudes, corrected for PSF losses, 
obtained from PSF-matched images in all but the Spitzer bands. 



The AGN and Galaxy Evolution Survey (AGES; iKochanek et al.l 120121 ) obt ained deep 
optic al spectra of approximately 25,000 sources in the Bootes field with Hectospec ( jFabricant et al 
20051 ) at the Multiple Mirror Telescope (MMT). The survey is designed to be statistically 
complete for several different samples limited to / < 20 for galaxy candidates and I < 22.5 
for AGN candidates. AGES is highly complete for AGN candidates to I < 21.5. The AGN 
candidates were targeted by their X-ray, radio and mid-IR properties, b ut not by t heir optical 
colors. This ensures that none of the optical selection biases (see, e.g., lFanlll999[ ) are propa- 
gated into the sample. We complement the AGES spectroscopy with ~ 2000 deeper optical 
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spectra from various sources, primarily from Keck (e.g.. lEisenhardt et al.ll2008l ). Since these 
data do not have a uniform selection function, they will only be of limited use in our analysis. 



2.2. WISE Observations 

The WISE mission observed the full sky in four mid-IR photometric bands with FWHM 
of 6" in Wl-3 and 12" in W4. We use the WISE all-sky data release, which includes all 
observations obtained during the fully cryogenic mission. WISE surveyed the sky in a polar 
orbit with respect to the ecliptic, simultaneously obtaining images in all four bands. Hence, 
the number of observations in a field increases with its ecliptic latitude. While fields near 
the ecliptic were typically o bserved 12 times, t he number can grow to several hundreds 



near the ecliptic poles (e.g., iJarrett et al.ll201ll ). The median coverage across the sky is 



approximately 15 frames per passband. In particular, the COSMOS field was observed with 
a median coverage of 11 frames per pass band, well below the median sky coverage. Detailed 



accounts of the mission are presented by IWright et al.l (120101 ) and in the WISE all-sky data 
release explanatory supplement^. 

The NDWFS Bootes field is at an ecliptic latitude of 46 deg, and hence WISE obtained 
an average coverage of 30 frames in each band, reaching 10cx depths in Wl, W2, W3 and W4 
of approximately 17.12 mag, 15.73 mag, 11.55 mag and 7.83 mag. For all sources, we use 
fluxes obtained through profile fitting. We limit the sample to S/N > 3 in Wl and W2, or 
equivalently to Wl<18.50 and W2<17.11. We match to other sources in Bootes by finding 
the closest IRAC [4.5] source within 2" with the constraint that no WISE (IRAC) source 
is matched to more than one IRAC (WISE) source. This results in a sample of 111,720 
matched sources. We note that the WISE magnitude limit is applied after cross-matching 
with the IRAC sources. 

Detailed comparison between WISE and Spitzer IRAC photometry has shown that 
WISE profile-fitting fluxes are typically underestimated for faint sources, and that the mag- 
nitude of the effect increases with decreasing IRAC flux, reaching offsets of a few tenths 
of a magnitude for the fainter sources (see section VI. 3 of the WISE all-sky data release 
explanatory supplement for details). While there is no simple prescription to mitigate it, 
this bias is unlikely to affect the results of our SED fits, as it is only significant for faint 
sources for which the deeper IRAC SDWFS magnitudes dominate the \ 2 °f the fit in the 
mid-IR. Because of this, we do not attempt to compensate for this WISE calibration issue 
in our study. 



i- 



http : //wise2 . ipac . caltech.edu/docs/release/allsky/ expsup/ 
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2.3. Spectral Energy Distribution Modeling 



We rely on SED modeling both to obtain physical insight into our AGN candidates, 
and to obtain photometric redshift (z p hot) estimates for all objects without available spectro- 
scopic redshifts. To fit the SEDs we use the non-n egative basis of low -resolution, UV through 
mid-IR SED templates for AGN and galaxies of lAssef et al.l (120 lOh . The basis consists of 
four empirically derived SED templates, where every object is modeled as a non-negative 
combination of the three galaxy SED templates (roughly corresponding to E, Sbc and Im 
types) and the single AGN template. For the AGN template alone, we allow reddening with 
a strength parametrized by E(B — V). For high-r e dshift sources, we model the intergalac- 
tic medium (IGM) absorpti on following iFan et al.l ( 120061 ) for Lya and Ly/3 absorption and 
Stengler-Larrea et al.l (ll995f ) for Lyman limit systems. The strength of the IGM absorption 
can also be fit beyond the standard mean absorption law, although this extra degree of free- 
dom often has a negative impact on the accuracy of photometric redshifts. A weak prior is 
used to keep E(B — V) as small as possible with the secondary effect t hat obscurat i on va lues 
may be slightly underestimated in some cases. We refer the reader to lAssef et al.l (120 lOf ) for 
further details. 



We follow the prescription detailed in lAssef et al.l (120101 ) to obtain photometric redshifts 
and fit the SEDs. Since photometric redshifts using 2 4/xm photometry have lower accura- 
cies when using these templates (see lAssef et al.ll2010l . for details), we derive photometric 
redshifts for all objects in the sample using only the broad-band photometry from FUV to 
W3. Adding MIPS 24/zm and W4 photometry, however, does not qualitatively alter our 
results. We assume the standard mean IGM absorption and use a luminosity prior for the 
galaxy co mponents based on the Las Campanas Redshift Survey (LCRS) r— band luminosity 
function ( ILin et al.lll996l ). We discuss the precision of the photometric redshifts in §4.11 



After obtaining photometric redshift estimates, we re-fit the SEDs of all objects now 
including the W4 and MIPS 24/im channels, and also fit for the strength of the IGM absorp- 
tion. Whenever possible, we use spectroscopic redshifts (z spec ). This approach ensures that 
we get the best SED model possible for each object. Several authors have determined that 
photometric redshifts for Type 1 AGN obtained solely with broad-band filter photometry can 
be w ildly inaccurate (see, e.g., iRowan- Robinson et al.ll2008l ; ISalvato et al.l 12009c lAssef et al. 
20101 ); this is discussed further in §4.11 in the context of our study. However, our spectro- 
scopic data is particularly deep and complete for AGN (see §2.ip . somewhat mitigating this 



issue. 



In order to reliably separate AGN from inactive galaxies, we use the parameter 

* _ £agn 



^host 



+ L 



AGN 
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where the luminosities correspond to the integrated specific luminosities of the best-fit tem- 
plates over the 0.1 and 30/zm wav elength range for the AGN template and 0.03 to 30 /zm 
for the host galaxy templates (see lAssef et al.l I2010L for details). The specific luminosities 
are calculated after correcting the AGN component for re ddening. W e refer to these as 
bolometric luminosities for the rest of the paper. Note that lAssef et al.l ( 120101 ) determined 
that a is insensitive to photometric redshift uncertainties as long as enough data exists to 
constrain the fit, in the sense that a can still be accurately determined for objects where 
reliable photometric redshifts cannot be measured. That is to say, while it is challenging to 
measure accurate photometric redshifts for AGN, particularly for Type 1 AGN, we are able 
to accurately disentangle the relative fractions of starlight and nuclear emission even when 
the redshift estimate is significantly in error. 



WISE AGN Color Selection 



In this section we stud y the com pleteness and reliability of WISE AGN selection. First, 
we discuss the criterion of iPaper t applied to the deeper WISE data in the Bootes field, 
while in §3.21 we improve our method by also considering the observed W2 magnitude of the 
sources. In §3.31 we compare our new method with others in the literature. 



3.1. Magnitude-Independent AGN Color Selection 



In|Pap_er_J we investigated the distribution of quasars in WISE color space in the COS- 
MOS survey field. We found that for objects with W2<15.05 mag (W2 S/N > 10 at that 
ecliptic latitude), the simple color cut based on the two shortest wavelength WISE bands 



Wl - W2 > 0.8, 



(2) 



provides an effective criterion to sep arate AGN f r om in acti ve galaxies. When compared to 
the IRAC color selection method of IStern et al.l (120051 . see lAssef et al.l |2010| for discussion 
about its reliability), the WISE color criterion selection recovers 78% of the IRAC-selected 
AGN with a 95% reliability. Notably, six of the AGN candidates selected by IRAC and 
WISE were not d etected in the 1.8Ms Chandra survey of the C OSMOS fi eld (C-COSMOS; 
Elvis et al.l 120091 ). suggesting they may be Compton-thick (see IPaper J for details). The 
reason behind the success of this criterion is clearly illustrated in Figure [TJ This figure shows 



mid-IR color as a function of redshift for the AGN and galaxy SED templates of lAssef et al. 
(120101 ). For z < 3, the W1-W2 color of unobscured AGN is well above the color cut of 



0.8 mag. At higher redshift, reddened AGN (E(B — V) > 0.4) can also be redder than this 
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color cut. In practice, however, it is exceedingly uncommon to find high redshift, highly 
reddened AGN bright enough to be detected by WISE and so much m ore lumi nous that 
their host galaxy to dominate the rest-frame optical emission. Indeed, in iPaper II we found 
no quasars with z > 3 in the sample. We also find very few galaxies with z > 1, since the 
WISE observations of the COSMOS field are not deep enough to find normal galaxies at high 
redshifts. While other populations such as some ULIRGs and brown-dwarfs can have even 
redder W1-W2 colors, they are too rare in comparison to AGN to be a significant source of 
contamination. Note that Figured] implies that this W1-W2 color selection is biased against 
AGN which are faint with respect to their host galaxies. If the flux in the WISE bands is 
dominated by the galaxy, the colors will drop below the selection limit, moving towards the 
galaxy locus at Wl — W2 ~ 0. We discuss this further in §4.31 



The criteria of IPaper II are readily applicable to the all-sky WISE survey, and are demon- 
strated to be both reliable and complete to the shallow depth of the WISE observations of 
the COSMOS field. However, due to the limited size of that field, it does not have the sta- 
tistical power to address many interesting and pressing issues in AGN studies, such as AGN 
evolution, accretion rates, and dust distributions. More importantly, since most of the WISE 
survey area has deeper coverage than in the COSMOS field, alternative selection criteria are 
valuable for a census of WISE-selected AGN in these deeper regions. IJarrett et al.l (120111 ) 
has shown that in the deepest WISE fields at the ecliptic poles, where Wl and W2 are 
confusion limited, the addition of W3 is a very useful aid in the identification of AGN. Our 
intention is to bridge these two extremes, proposing a robust WISE AGN selection technique 
for fields with intermediate depth. 



In order to extend the study of IPaper II . we turn to the NDWFS Bootes field, which 
helps with both issues highlighted above: it has a WISE median coverage of 30 frames, 
almost three times that of COSMOS, and ex tends oy er 9 deg 2 , an area 4.5 times larger. We 
start by replicating the selection criterion of IPaper II in the Bootes field, but up to the W2 
10a depth provided by the full co-added data. While this implies a sample 0.68 mag deeper, 
it maintains the error properties of the sample and so provides a meaningful comparison. 

The left panel of Figure [2] shows the WISE color distributi on of sources in ou r sample. 
Comparing to sources selected as AGN by their IRAC colors (jStern et al.l 120051 ). we find 
that a simple Wl-W2>0.8 color cut identifies 70% of the IRAC-selectecl AGN with 70% 
re liability. Compared to applying this cut at the 10cx WISE depth of COSMOS as reported 
Paper II . the drop in completeness is relatively small, from 78% to 70%. The decrease in 



m 



reliability from 95% to 70% is very significant, however, and is simply due to the modest 
increase in field depth. If we limit the Bootes field analysis to the 10a W2 level of the 
COSMOS field (W2<15.05, right panel of Fig. [2]), we recover similar AGN demographics to 
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that reported in iPaper J . with 78% completeness and 94% reliability. 



The lower completeness means an increase in IRAC-selected AGN detected but not 
identified by the simple WISE color criterion. This is likely due to a combi nation of (i) a 
small number of z > 3 Type 1 AGN, which are known to be excluded by the IPaper t color 
selection; (ii) low-redshift, low-luminosity AGN with hosts bright enough to mov e their mid 



IR co lors below the WISE selection limit, but red enough to be picked by the IStern et al. 
( 120051 ) IRAC selection criteria; and (iii) a higher incidence of contamination by z ~ 0.5 star 



forming galaxies to the IRAC selection criterion which artificially lowers the completeness 
- t hough we note that this contaminatio n is expected to be small at the depth of SDWFS 



see 



Assef et al.ll2010l : iDonlev et all 120121 ) 



The cause of the significantly lower reliability obtained in the Bootes field compared 
to the COSMOS field is readily apparent in Figure [2J The modestly deeper WISE sample 
increases the number of contaminating galaxies, particularly to the left of the QSO locus. 
These correspond to high-redshift (z ~ 1 — 1.5) galaxies. The observed W2 magnitude of 
z ~ 1-1.5 L» galaxies evolves very slowly with redshift (see, e.g., Fig. 1 of lEisenhardt et al. 



20081 ). so there is a huge increase in contaminatio n as soon as the W l magnitude limit is deep 
enough to begin including these galaxies. In the IStern et al. (b()05l ) IRAC selection criterion 
this problem is controlled using the [5.8]-[8.0] color, but the longer wavelength WISE bands 
are too shallow to help. 



3.2. Magnitude-Dependent AGN Color-Selection 



It is apparent from Figure [2] that an improved method to select AGN may be possible if 
we allow our color cut to evolve with magnitude since the major contaminants are either low- 
redshift, nearby star-forming galaxies which are intrinsically faint, or high-redshift, passive 
galaxies that are luminous enough to be bright in the WISE bands. In order to design a 
magnitude-dependent AGN color selection method that is applicable over the whole sky, we 
will go to fainter WISE fluxes than afforded by the 10cx W2 limit in the Bootes field. Note, 
however, that as we go to fainter W2 magnitude s, it becomes unre liable to use an AGN 
control sample based on the IRAC color criteria of IStern et al.l (120051 ) . as it is susceptible to 
contamination by high-redshift g alaxies once the err ors in the SDWFS IRAC [5.8] and [8.0] 
fluxes become too large (see, e.g.. lDonley et al.ll2012l ). Instead, we define the control sample 
as all objects whose best-fit UV-mid-IR SEDs have a strong AGN component, as indicated 
by requiring a > 0.5. For significantly lower levels of AGN activity, it becomes necessary to 
differentiate between objects where the AGN component of the fit is real and when it has 
only been used to mathematically improve the \ 2 t° accommodate lower quality photometry 
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or mimic a galaxy component missing from the templates. This f alls beyond the scope o f 



the current work, and a full analysis on this topi c is pr esented by IChung et al.l ( lin prep.l ). 



Assef et al.l (120101 ) has shown that the lStern et al.l (120051 ) criterion is biased towards objects 



with large a values (see also §4.3p . so we are not considerably changing the physical properties 
of the control sample by using this definition. 

Figure [3] shows the completeness and reliability obtained as a function of W2 magnitude 
and the minimum W1-W2 color limit adopted to select AGN. We have required a minimum 
detection threshold of 3a for Wl in order to have a reasonably precise WISE color. At 
bright W2, a color cut of 0.6 is sufficient to obtain high reliability and high completeness. 
Towards fainter W2 magnitudes, high reliability requires redder color cuts in order to remove 
contaminating galaxies, which also leads to lower completeness. The completeness of a color 
cut is relatively independent of W2 magnitude. 

Figure [3] shows the bluest W1-W2 color at which 90% and 75% reliability is reached for 
a given W2 magnitude. While there is significant noise in these curves, they are reasonably 
well described by an exponential in W2 2 . Hence, we propose a WISE AGN color selection 
limit optimized for reliability given by 

Wl - W2 > a R exp {/3 R (W2 - 7r ) 2 } . (3) 

For W2<17.11, we achieve a reliability of ~90% with (a mo , (3 R90 , 7rao ) = (0.662,0.232,13.97). 
The corresponding values for a reliability of ~75% are (ai?75, /3r75, 7^75) = (0.530, 0.183, 13.76). 
The 90% (75%) reliability criterion reach our imposed Wl S/N > 3 limit at a W2 magnitude 
of 16.26 (16.45). Only considering objects brighter than this limit in W2, the 90% reliability 
criterion identifies 1174 AGN candidates, of which 1060 (90%) have their bolometric lumi- 
nosities dominated by the AGN (e.g., a > 0.5). For the 75% reliability criterion, we identify 
2306 AGN candidates, of which 1752 (76%) are AGN-dominated. At a shallower depth of 
W2 < 15.73, corresponding to S/N > 10 in the Bootes field, the 90% reliability curve identi- 
fies 1051 AGN candidates of which 950 (90%) are AGN-dominated, while the 75% reliability 
curve identifies 1582 AGN candidates of which 1200 (76%) are A GN-domi nated. For com- 



parison, a simple Wl-W2>0.8 cut (e.g., similar to the criterion of iPaper 11 . but without its 
magnitude cut) finds 1746 AGN candidates with 74% reliability at a depth of W2 = 15.73. 
This census increases to 17997 AGN candidates to a depth of W2 = 17.11, albeit with a 
reliability that drops to 45%. These statistics are summarized in Tables [T] and [2J In the next 
section ( §3.31) we compare this new magnitude-dependent AGN selection criterion to several 
other WISE AGN selection criteria that have recently been proposed in the literature. 

Figure [3] also shows the reddest W1-W2 color at which 90% and 75% completeness 
is reached for objects with bolometric luminosities dominated by the AGN emission. In 
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this case, a reasonable description of the completeness boundary is given by the magnitude 
independent color cut 

Wl - W2 > 6 C , (4) 
where <5c90 = 0.50 for 90% completeness and 6c7§ = 0.77 for 75% complete ness. No te that 



the 75% completeness criterion is basically equivalent to the cut proposed by iPaper II . shown 
in equation [2], but without the flux cut. It is important to stress that these criteria are 
appropriate only for strong AGN with respect to their hosts due to the a > 0.5 requirement. 

While the magnitude dependence in the reliability-optimized criterion is caused in part 
by the much higher number of high redshift galaxies at fainter fluxes (see discussion in §3.ip . 
it is also driven by the increasingly large errors in Wl and W2 at lower S/N. Since the 
S/N of WISE observations varies significantly across the sky, in principle ccr, (3r and 7r 
may depend on ecliptic latitude. To test the strength of this dependence, we simulate the 
distribution of W1-W2 and W2 magnitudes for different WISE field depths and estimate or, 
/3r and 7r in each case. For this we use the magnitudes obtained from the SED modeling 
of every object and we approximate that the S/N for a flux F depends on field depth as 



S „ N F F 

where Np is the number of individual WISE lis frames used to build the catalog image, K% is 
a constant and F^ is the background flux. This formulation neglects the effects of confusion 
as well as the variation of F^ Y with sky position is not uniform across the sky. However, 
our approximation should give a good general idea of how the parameters in question vary 
with N F . As discussed earlier, jVj? ootes = 30, and we use the Bootes data to fit for K\ and 
-Psky We find that from Np = 10 to 50, no significant variation is observed for the or and 
7r parameters for both the 90% and 75% criteria. The parameter (3r is observed to decrease 
linearly by a factor of ~ 4 between Np = 8 and Np = 25, and is approximately independent 
of the depth of the field for N F > 25. Considering then a modified 0' R = (5.41-0.176^) 
may be necessary to achieve the proposed reliability levels in WISE fields with Np < 25, 
where (3®°° is the value obtained for the Bootes field. We also repeat the experiment for the 
completeness optimized criteria, and find that in the same range of Np, 8c is approximately 
constant. 

We refer to the reliability-optimized selection as -R90 and -R75 for 90% and 75% relia- 
bility, respectively, while C90 and C75 refer to the completeness-optimized criteria. Whether 
scientific interest lies in maximizing completeness or reliability depends on the problem at 
hand. However, it will be most common to wish to maximize reliability, so in the next 
sections we will focus on results for the highest reliability selection. 
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3.3. Comparison with the Literature 



AGN identification using mid-IR broad-band photometry is now a well studied problem 
The first classification schemes were developed for Spitzer IRAC and MIPS photometry (see 



e.g. , lLacv et al.ll2004 120071 : IStern et al.ll200a lAlonso-Herrero et al.ll2006l : Messias et al.ll2012[ ) 
and have been shown to be very successful in terms of both reliability and completeness. 
In the previous section we developed four AGN identification schemes using WISE Wl and 
W2 photometry, optimized to produce samples with different levels of either reliability or 
completeness. Several other WISE criteria have also recently been developed. Here we briefly 
discuss several of these criteria and discuss how they compare to our selection criteria. This 
is meant to be an illustrative rather than an exhaustive exercise; we do not discuss all the 
published mid-IR selection techniques. 

Tables [fl [2] and [3] show the surface density of AGN candidates, the i r relia b ility and their 
compl etene ss for the selectio n criteria of 



Paper I. 




Jarrett et al. 


(2011 


)• 


Mateos et al. 



(|2012j) and lWu et al.l fl2012bf ). Table □ restricts the samples to W2<17.11 (W2 S/N > 3 in 
Bootes), while Table [2] uses the more restrictive flux cut of W2<15.73 (W2 S/N > 10 in 
Bootes). The samples used in all tables are also restricted to Wl<18.50, but no restriction 
is applied in W3 and W4. For completeness, Table [3] further limits the sample to W2<15.05 
(W2 S/N > 10 in COSMOS), which is representative of the shallowest WISE observations, 
but we will not discuss it in detail. We also include the selection criteria proposed by 



Assef et al.l ( 120101 ) . which was obtained by simulating WISE photometry using SED models 
of all objects in SD WFS. For comparison, we also show the numbers for the IRAC-based 



selection criteria of lLacv et al.l ( 12004 |2007J), IStern et al.l ( 120051 ) and Messias et al.l (120121 ). 
calculated using the SDWFS photometry. As was done in the previous section, reliability 
and completeness are measured against the number of objects whose SED fits have d > 0.5. 

For the WISE selection methods, regardless of the W2 depth, the most re liable sample 
is that based on the Wl, W2, W3 and W4 selection criteria of Assef et al.l (120101 ) . with a 97% 
and 98% reliability for W2<17.11 and W2<15.73, respectively. However, because it requires 
that W4 is detected, it also has the lowest completeness (9% and 19% for W2<17.11 and 
W2<15.73, respectively) as measured by the surface density of AGN candidates, with only 
43 deg -2 candidates with W2<17.11. Our R$ is the second most reliable criterion, with 
90% reliability by design, but it has a much higher completeness, with 53% for W2<15.73 
and 9% for W2< 17.11, which translates into AGN c andidate surfaces d ensities of 117 
and 130 deg -2 respectively . The four-band criterion of iMateos et al.l ( 120121 ) also has high 
reliability, although it is below our Rgo criterio n in both reliability a nd c ompleteness. Both o f 
the Wl, W2 and W3 based selection criteria of ljarrett et al.l ( 1201 ll ) and lMateos et al.l (120121 ) 
are similar in reliability and completeness, comparable to our Rj 5 criterion for W2<15.73, 
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but somewhat less reliable for W2<17.11. As discussed by lJarrett et al.l (120111 ). the strength 
of these criteria are in the deepest WISE fields, where Wl and W2 depths are below the 
confusion limit. 

With respect to th e IRAC -based criteria, the most reliabl e of t hose shown are the 
criteria of iMessias et al.l (120121 ). followed by that of IStern et al.l (120051 ) . In terms of com- 
pl eteness levels, we not e that all criteria shown are similar, except for the "KIM" criterion 
of IMessias et al.l (120121 ) based on K s , [4.5], [8.0] and MIPS 24/im photometry, which has a 
significantly lower completeness of 22%. It is also important to notice that, in principle, 
completeness and reliability could be improved by further adding more information based on 
other wavelength r e gimes . Such is the case, for example, with the "Six" selection scheme of 
Edelson fc Malkanl (boi2h . which combines WISE, 2MASS and ROSAT data to identify the 
brightest AGN in sky. Including this kind of selection is, however, beyond the scope of this 
comparison. 



4. Properties of WISE AGN Candidates 

In this section we study the properties of the WISE AGN candidates selected using 
the criteria developed in the previous section. We first discuss the accuracy to which we 
can determine photometric redshifts for them. In §4.21 we discuss their redshift distribution, 
and in §4.31 we discuss the parameters derived from our SED fitting. In §4.41 we present 
spectroscopic observations of a sample of photometrically selected high redshift Type 2 AGN 
candidates. 



4.1. Photometric Redshift Accuracy for WISE AGN Candidates 



Several authors (e .g., iBrodwin et al.l 120061 : iRowan-Robinson et al.l 120081 : ISalvato et al. 



20091 ; lAssef et al.ll2010l ) have shown that photometric redshifts of Type 1 AGN are relatively 
inaccurate when relying solely on broad-band photometry, as is our case. This is mostly 
due to the lack of strong spectral features that are necessary for anchoring the photometric 
redshift estimates. Our AGN sample is, however, brighter than those typically studied for 
photometric redshifts, and has a considerable number of Type 2 AGN. Photometric redshifts 
for Type 2 AGN may be better because the spectral features of the host galaxy are relatively 
stronger. 

We estimate photometric redshifts as discussed in §2.31 using, in addition to WISE, 
all the UV through mid-IR broad-band photometry of the field described in §2.11 As in 
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Assef et all (120101 ). we quantify the photometric redshift accuracy using the statistic 

1/2 



Az 



1 

N 



'phot 



spec 



1 + Z, 



spec 



(6) 



where the index i sums over all objects in a sample and N is their total number. This 
estimate of the dispersion, however, is typically driven by outliers, so we also estimate Azg§, 
the dispersion calculated including only the 95% of objects with the best photometric redshift 
estimates. 

Panel a) of Figure H] shows the spectroscopic and photometric redshifts obtained for the 
full W2 depth Rg AGN candidates, limited to objects with a > 0.5 to be certain we only 
study the objects of interest. Table @] shows the dispersion as well as the median offsets 
for the remaining criteria, again limited to a > 0.5. It also shows the number of AGN 
used to compute the statistic and the fraction of objects in every selection criteria that have 
spectroscopic redshifts. Irrespective of the selection method, the photometric redshifts are 
fairly inaccurate, wit h Az 9 s = . 20 — .23 (Az = 0.27 — 0.31). This is consistent with the 
results presented by lAssef et all (120101 ) for a similar, but fainter, sample of objects. The 



pile-up of objects at very low z phot is a degeneracy caused by the galaxy luminosity prior. 
However, these are only a small part of the sample, and eliminating the prior results in even 
less accurate estimates for the general population. Panel b) shows that little is gained in 
terms of the accuracy when limiting the sample to the brighter W2 S/N > 10 objects. The 
same is observed when the sample is further limited by requiring / < 20, as shown in Panel 



If we further split the bright, final sample (W2 < 15.73 and / < 20) and only investigate 
objects with considerable obscuration, E(B — V) > 0.5, we find that Az$ 5 drops by ~ 45%, 
although Az either decreases only slightly (Rqq and C75), or increase s {Rm and Can). T hat 
these accuracies are still much worse than the Azg 5 ~ 0.04 found by lAssef et all ( 120101 ) for 
galaxies of equivalent brightness is most likely due to the sample requirement that a > 0.5, 
meaning that even though reddened, the accretion disk emission is still likely dominant, or 
at least significant, in many of the broad-bands used. 



As mentioned earlier, lAssef et al.l (120101 ) showed that although photometric redshifts 
for AGN based on broad-bands can be inaccurate, the value of a obtained from the corre- 
sponding SED fit is insensitive to the redshift accuracy, i.e., that a is relatively independent 
of photometric redshift. As one of our goals is to study the obscuration in AGN in a sta- 
tistically significant manner, we can ask if this holds for the inferred reddening of the AGN 
component. Unfortunately this is not the case; while objects with good photometric redshift 
estimates have consistent estimates of E(B — V), objects where z p hot is considerably in error 
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have inaccurate and typically underestimated values of E(B — V). This will be of particular 
importance in §|5j 



4.2. Redshift Distribution of WISE AGN Candidates 

Using the cuts developed in the previous section we now study the redshift distribution 
of the different samples of AGN candidates. Although we have a large amount of spec- 
troscopic observations in the Bootes field, we are still missing spectroscopic redshifts for a 
considerable number of our AGN candidates (see TableHJfor details). For the objects without 
spectroscopic redshifts, we use the photometric redshift estimates detailed in §2.3[ although 
these may not be very accurate (see §4. ID . We focus on the Rgo sample, which mitigates this 
issue as these objects are the ones most likely to have spectra from the AGES survey (see 



2.11 and iKochanek et al.ll2012l . for details). 



Figure |5] shows the redshift distribution of the W2, Wl S/N > 3 depth R$ sample. 
The resulting distribution of AGN is double peaked, with the main peak at 1 < z < 2 and 
a smaller peak at z ~ 0.25. Almost no objects are at z > 3. This distribution reflects 
that WISE has a high sensitivity to obscured AGN at lower redshifts, where the AGN 
emission still dominates the observed Wl and W2 fluxes. However, as the redshift increases 
or the galaxy host contributions become larger, the bias against obscured sources increases. 
This causes the minimum at z ~ 0.75, followed by an increase simply from the increase 
in comoving volume probed. The W1-W2 color of Type 1 AGN is reddest at 1 < z < 2 
(see Fig. CD), and progressively gets bluer at higher redshift, falling completely out of the 
selection criteria by z ~ 3. A similar behavior is observed for the .R75, C75 and C90 samples, 
although since contaminants appear preferentially at high redshifts (see §3.2p . the balance 
between the peaks for the complete a samples is modified. Figure [5] shows that a large 
number of the -R90 AGN candidates (56%) have spectroscopic redshifts. Furthermore, objects 
lacking spectroscopic redshifts tend to follow a similar redshift distribution, implying that 
although the uncertainties in the photometric redshifts are very large, they do not seem to 
systematically bias the distribution. Limiting the samples to only the brighter W2 S/N > 10 
objects does not significantly change the shape of the redshift distribution. 



4.3. SED Analysis of WISE AGN Candidates 



A simple way of quantifying the contamination rates in the criteria we have defined is by 
looking at the best-fit combination of SED templates to their photometry. The most relevant 



- 16 - 



parameter is a, defined in equation ([!]). As mentioned earlier, this parameter ha s the useful 



prop erty of being relatively insensitive to photometric redshift uncertainties (see lAssef et al. 
201(1 for details). Figure El shows the distribution of a for our full-depth R 90 and C 90 AGN 
candidate samples. We find that the R 90 sample is skewed towards objects dominated by 
their AGN component, with almost no objects being best-fit as inactive galaxies. This 
feature is also observed, although to a somewhat lesser degree, in the R75 sample. The C90 
sample, on the other hand, shows a very considerable peak at a = 0, as expected given its 
low reliability but high completeness. It also shows, however, a very significant increase in 
the number of objects with intermediate a values. These dominate the distribution for a > 0. 
Most such objects probably correspond to real AGN with high host fractions, implying that 
our reliability optimized criteria is strongly biased against such objects. 

It is well known that the luminosity of the spheroidal component of the host galaxy is 
correlated with the mass of its central SMBH, and that this relation is roughly linear: Lhost ~ 
Mbh (see, e.g.. iMagorrian et al.l 1 19981 ; iFerrarese fc Fordl 120051 ; I Graham! 120071 . although also 



see iGrahamll2012l for possible deviations). Some authors have postulated that the correlation 
is also present, and non-evolving, w hen considering the total host galaxy luminosity instead 
of just the spheroidal component (IBennert et al.ll2010l ). Regardless, since the Eddington 
luminosity LEdd is directly proportional to Mbh, the Eddington ratio ^Edd can be expressed 
as 

„ Lagn L agn a , . 

«Edd = - T ~ 7 = ~ r- (7) 

-^Edd -^Host 1 — a 

Hence, to first order, AGN whose bolometric output is dominated by the AGN emission 
(i.e., have high a values) also correspond to objects emitting at a high ^Edd- Similarly, those 
galaxies for which stellar light represents a higher fraction of their total bolometric output 
(i.e. low a) are likely radiating at lower Eddington ratios. So, in a physical context, we see 
that our reliability-optimized selection criterion is strongly biased against AGN radiating at 
low ^Edd, but as our selection criterion is shifted to emphasize high completeness, we start 
recovering them. 

Figure E] also shows the distribution of a limited to W2 < 15.73. While the R 90 sample 
looks nearly the same, the C90 sample exhibits a different distribution, with the contamina- 
tion (a = peak) and the skewness of the distribution shifting to be more similar to the Rgo 
sample. Partly, this is because the largest contamination in completeness-optimized samples 
comes from high-redshift galaxies, which are avoided by the R$o criterion. That the peak at 
high a is increased is possibly due to low Eddington ratio AGN simply being less luminous 
on average. 

The other SED-fit parameter of interest is the amount of obscuration towards the AGN, 
which is shown in Figure [7] for the Rqo AGN sample at full and 10cr W2 depths. The most 
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important result to notice is that the WISE AGN selection is sensitive even to objects with 
high obscuration. In order to interpret the distribution, however, we need to deal with two 
issues. The first is that because of the algorithm design, the reddening may be slightly 
underestimated. Second, we need to take into account the selection function of AGES, 
since the reddening obtained from objects with only photometric redshifts estimates can be 
inaccurate ( §4.ip . We deal with both issues and present a detailed study of the reddening 
properties of AGN in Section [51 



4.4. Keck Observations 

To highlight the power of WISE in finding highly obscured quasars, we obtained ad- 
ditional spectroscopy of 12 AGN candidates at the Keck Observatory in April 2011. Since 
the AGES spectroscopy is limited to / < 22.5, and is highly complete for / < 21.5, we 
emphasized opticall y fainter candidates which are bright in W2, selected on the basis of an 



early version of the iPaper 11 criteria. We furthermore required them to not have a measured 
redshift. Because these observations used a preliminary version of the WISE data, some of 
the WISE colors and positions changed relative to the more accurate all-sky release. There- 
fore, the sources we observed have a range of WISE colors and optical magnitudes; based 
on the WISE all-sky data release, all but one of them are sufficiently red in Wl— W2 to be 
classified as AGN candidates by at least the Cqq criterion, but not all of them are optically 
faint (/ > 21). However, most of the targets do meet the Rgo selection criterion and prove 
to be bona fide obscured quasars (see Table [5] for W1-W2 color, W2 magnitude and AGN 
classification criteria met by each target). 

We observed three Keck slit masks in the Boot es field with the DEep Imaging Multi- 



Object Spectrograph (DEIMOS; iFaber et all 120031 ) on UT 2011 April 1-3. We used the 



4000 A order-blocking filter and the 600 £ mm 1 grating (blazed at 7500 A; resolving power 
R = A/ A A ~ 1600 for the l'/2 wide slitlets we employed). We observed a single additiona l 



mask using the dual-beam Low Resolution Imaging Spectrometer (LRIS; lOke et al.lll995f ) 
on UT 2011 April 28. The LRIS observations employed the 300 I mm^ 1 grism on the blue 
arm of the spectrograph (blazed at 5000 A; R ~ 500), the 400 i mm -1 grating on the red 
arm of the spectrograph (blazed at 8500 A; R ~ 700), and 6800 A dichroic. Data reduction 
followed standard p r ocedu res, and we flux calibrated the data using standard stars from 



Massev & Gronwalll fll990h . 



Table [5] summarizes the results for these observations, including measured redshifts, 
selection criteria, and the best-fit a and AGN reddening parameters for the adopted redshift, 
with errors derived from Monte Carlo re-sampling of the data. The Appendix presents 



-18 - 



the results for additional Bootes targets observed on these masks. We include the quality 
("Q") of each spectroscopic redshift. Quality flag "A" signifies an unambiguous redshift 
determination, typically relying upon multiple emission or absorption features. Quality flag 
"B" signifies a less certain redshift determination, such as the robust dete ction of an isolate d 
emission line, but where the identification of the line is uncertain (e.g., IStern et all 120001 ) . 
Quality flag "B" might also be assigned to a source with a robust redshift identification, but 
where some uncertainty remains as to the astrometric identity of that spectroscopic source. 
We consider the quality "B" results likely to be correct, but additional spectroscopy would 
be beneficial. 

Figure [8] shows the best-fitted SEDs for each of the 8 targets, from the original 12, 
whose all-sky release WISE Wl— W2 colors classify them as AGN by either the R 90 or R 75 
criterion. Upon inspecting the optical images, we believe the bright, discrepant J-band flux 
of W1430+3530 is most likely due to a bright star within 30" contaminating the photom- 
etry. Although most of these objects appear to be real AGN based on their broad-band 
SEDs, many lack strong, high-ionization lines such as C IV, Mg II and Ne v, even though 
lower ionization lines common for star-formation are indeed observed (see Table The 
X-ray communit y has noticed a relat ed population of X-ray bright, optically normal galaxies 
(XBONGs; e.g., ICivano et all 120071 ) where the X-ray luminosities require the presence of 



an actively accreting SMBH while optical spectroscopy reveals an apparently normal, inac- 
tive galaxy. Several explanations have been offered to explain such sources, ranging from 
systematic effects that dil ute the AGN signa ture for the wide slit widths typically used for 



these distant sourc es (e.g.. iMoran et al.ll2002[ ). to radiatively inefficient accretion flows (e.g. 



Trump et eil.l l201lh . Alternatively, at least some of these objects could be better described 
as AGN dominated LIRGs or ULIRGs, where the lack of high ionization emission lines and 
the red host color may be explained by large scale obscuration. Some evidence of the Si 
9.7/zm absorption feature typical of ULIRGs may be present in a few cases (W1427+3408, 
W1431+3525, W1432+3523 and W1432+3 526), although thi s feature may also be observed 



in AGN under certain conditions (see, e.g.. lFeltre et al.ll2012l ) 



5. Dust Reddening in AGN 

In this section we study dust obscuration properties of a set of 362 z < 1 AGN well 
detected by WISE with spectroscopic redshifts, I < 20 and a > 0.5 in the Bootes field. As 
argued earlier, low redshift (z < 1) WISE AGN selection criteria are relatively insensitive to 
obscuration since they rely on the hot dust emission from the dust torus instead of on the 
blue colors of the unobscured accretion disk emission, as per optical selection. Hence, we 
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can use WISE to study the properties of dust obscuration in AGN. 
AGN unification models (see, e.g. 



Antonuccilll993t lUrry fc Padovanilll995l ) propose that 



Type 1 and Type 2 AGN are physically equal but are observed at different inclination angles 
relative to the obscuring material near the AGN. Typically it is assumed that the accretion 
disk, responsible for the A < 1/zm continuum emission, extends to radial scales of ~ 10 AU, 
and is surrounded by highly ionized gas responsible for the broad emission-lines. On larger 
scales (~ 1 pc) there is dust in the general shape of a "torus" or a flared disk, responsible 
for the Type 1/Type 2 dichotomy, that absorbs the optical radiation from the accretion 
disk and re-emits it in the mid-IR. We refer to this structure as the torus, as is commonly 
done, although we do not a priori assume a shape for it. The inner edge of the dust torus 
is determined by where the dust reaches its sublimation temperature due to heating from 
the accretion disk. Such hot dust produces the emission observed to dominate the mid-IR 
portion of the AGN SED. Furthermore, narrow emission line s are observed in both Type 1 
and Type 2 AGN, and are known not to be polarized (e.g., lAntonuccil Il993l ) . so the dust 
structure must be smaller than the narrow line region (~ 1 kpc). Given that the Type 
1/Type 2 dichotomy is also manifested in the neutral hydrogen absorption of the X-ray 
emission, the torus must also be associated with the absorbing gas. 

Many pr operties of the dust torus have been extensive l y studied. For example, several 



authors (e.g..lKrolik fc Begelmanlll988l : iNenkova et al.ll200l 120081 : Elitzur fc Shlosmanll2006 



Tristram et al.ll2007| ) have argued that the dust in the torus must be in optically and geomet- 



rical 



2005 



y thick clumps to reproduce observations, while others (e.g.. lDullemond &: van Bemmel 



Fritz et al. 



2006) argue the dust may be smoothly distributed. A recent study by 



Feltre et al. J2012h ~ however, suggests that given the same dust composition and the same 
illuminating source, the difference in the broad-band shape of the SEDs from these dust 
configurations may be too subtle to distinguish between both scenarios with current data. 
In a more global sense, the geometry and evolution of the obscuring structures have also been 



Ueda et al. 


2003). 


Simpson 


(2005) 



has shown using the Type 1 and Type 2 AGN optical luminosity function from SDSS that 
the fraction of Type 2 AGN increases with decr easing accretion disk luminosity , and a similar 



beha vior has been observed f or radio galaxies ( Lawrence 



200J) and in the X-rays (e.g.. lUeda et al.ll2003 



Hasinger 



199 



Simpson! Il998t iGrimes et al. 



20041 ) . Such a behavior can be nat- 



urally expected if the scale height h of the obscuring material is independent (or not linearly 
related) to the radial size of the structure (R oc a/Xagn), such that for brighter AGN the 
dust effectively covers a smaller solid angle as viewed from the SMBH. Th is scenario is usu - 
ally referred t o as the " reced ing torus model", and was first proposed by iLawrencd (Il99ll ). 
In particular, I Simpson! (120051 ) has shown that observations appear to be best reproduced if 
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a OC 1^ AGN - 

Combining WISE and all the ancillary observations in the NDWFS Bootes field and 
performing the SED modeling as detailed in §2.3} we can study the average properties of 
dust obscuration in AGN by counting the number of objects observed per unit reddening. 
This approach allows us to quantify the fraction of AGN that can be classified as Type 1, and 
also, in principle, to differentiate between different dust geometries and compositions. Note 
that obscuring dust may also be present in the interstellar medium (ISM) of the respective 
host galaxies, and while we expect AGN obscuration to be mainly driven by the dust in 
the torus and hence delineate the discussion in that direction, we further address galactic- 
scale obscuration in §5.31 Since our sample is, in essence, flux limited, our analysis must 
properly take into account its selection function. Specifically, it must account for all the 
biases against highly obscured objects, since higher AGN obscuration can also make the 
objects appear much fainter depending on the relative contribution and S EP shape of the 



host galaxy. Fortunately, the SED fitting approach of lAssef et al.l (1201 Ol see also §2.31) is 
well suited to assess and correct for our survey incompleteness. In the next section we detail 
our sample selection function. In §5.21 we detail the formalism we use to incorporate the 
selection function in our measurement of the reddening distribution, while in §5.31 we show 
and discuss the resulting distributions. 



5.1. Sample Selection Function 

To study the reddening distribution, we use a subsample of the larger sample described 
in §2J We require that objects have W2 < 15.73, a measured spectroscopic redshift such that 
E(B — V) is accurately estimated (see §4.ip . and a > 0.5 to minimize possible non-AGN 
contaminants. Note that we do the initial selection with the SED fits obtained including 
all the priors described in §2.3[ which is necessary to ensure all obscured AGN are real. As 
discussed there, this can lead to slightly underestimated AGN obscuration. Hence, once the 
sample is selected, we re-fit the SEDs removing all priors described in §5.31 to obtain the 
final E(B — V) values, although we note not removing the prior does not qualitatively affect 
our results. We visually inspected the SED fit of every source and eliminated 16 galaxies 
where we believed the AGN classification or the reddening values were spurious due to bad 
photometry. 

We further require that the redshift was determined by AGES, since its well determined 
selection function is a crucial component of our analysis. Since many of the objects we 
consider are extended in the NDWFS imaging and were not necessarily targeted as AGN 
candidates by the AGES survey, we must also restrict our sample to objects with / < 20, 
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resulting in a final sample size of 362 AGN. AGES was designed to ensure subsamples are 



statistically complete to / < 20 for galaxies and / < 21.5 for AGN (see lKochanek et al.ll2012 
for details). In order to do this, AGES used a sparse sampling algorithm for galaxies, such 
that for every defined galaxy subsample, a spectrum was attempted for all objects brighter 
than a certain magnitude limit and for a percentage (typically 20-30%) of randomly selected 
fainter galaxies down to a certain magnitude. For example, the main J-band selected galaxy 
sample was observed in full for / < 18.5 and 20% of the galaxies were followed in the range 
18.5 < I < 20. In contrast, there was no sparse sampling for AGN candidates as AGES 
attempted to get spectra of all of them. Every subsample was assigned a selection code, 
where -P™ sparse is the probability that object i of the subsample with selection code n was 
selected for spectroscopy due to the sparse sampling algorithm. 

In addition, the fraction of sources with a successfully measured redshift depends on 
/-band magnitude. While the survey design minimized the magnitude dependence beyond 
the sparse sampling, there is still a dependency simply because it is more difficult to obtain 
redshifts for fainter sources in a fixed integration time. Using the full results of the AGES 
survey, we estimate for every selection code the fraction of objects for which spectra were 
attempted and a redshift was measured as a function of /-band magnitude, P™ 7 . 

In order to correct for the selection function, we need to estimate for every object the 
probability that objects with the same optical and IR magnitudes would have been observed, 
so that we can statistically account for those without spectroscopic observations. Since every 
object may have been targeted for more than one of the different subsamples, we need to 
consider the joint probability of all subsamples. Let P™ = P™ sparse x P™ 7 and let N be 
the total number of subsamples object i is part of. We define C(N,k) to be the sum of 
all possible combinations of k element products of the P™ terms, such that, for example, 
C(3, 1) = Pi + Pf + Pf, C(3, 2) = PfPf + PfPf + PfPf and so on. The probability a 
spectroscopic redshift would have been obtained for objects like object i is then given by 

N 



£(-l) fc+1 C(N, k). 



k=i 



It can be shown that if any of the terms P™ = 1, then Pj = 1, as would be expected. For 
consistency with the original AGES selection, we use the original catalogs of AGES to assess 
the spectroscopic completeness rather than the catalogs described in $2j 
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5.2. Method 

In order to incorporat e the selection function, we adapt the step-wise maximum like- 



lihood method (SWML) of lEfstathiou et al.l (119881 ). For the remainder of this section, we 



define E BV = E(B — V). Our goal is to estimate the distribution ^(E BV ) = dn/dE BV , where 
we remind the reader that E B y corresponds to the reddening only over the AGN component, 
not over the host galaxy (see §2.31 for details). The probability of finding an object with a 
given reddening E BV is given by 



Pi oc I E . yM f (EW) ) . (9) 

Jo sv ' Max i(E BV )dE BV ^ 

where Cj = Pf 1 is calculated using equation (JSJ) and E BVMax is the maximum reddening 
object i could have and still be in our sample, which we detail below. We estimate £W,Max by 
varying the reddening of the AGN component of the best-fit combination of SED templates 
but keeping the amplitude of the components fixed (see §2.31 for details). 

To apply the SWML method we discretize the function £(E B y) in bins of E B y rather 
than assuming a parametric form. We divide ^{E B y) into N p bins of value centered at 
reddening values with widths AE BV . We can then rewrite equation ([9]) as 

p . oc ( ^W(E BV -e k ) \° l 

% KE^H^-E^J^AEzy) ' 

where 

w{x) = / 1 'f -AW2<*<AiW2, (n) 
I otherwise 

if x < -AE BV /2, 

HU) = { - if - AE BV /2 < x < AE BV /2, (12) 

if x > AE BV /2. 

The likelihood C of our sample being drawn from the distribution t^(E B y) corresponds to 
the multiplication of the values of all JVagn objects in our sample. Taking the gradient of 
C with respect to £ = (£i, ...,£&), we can find that the values that maximize the likelihood 
are given by 

p Ap ES GN Q W(E bv - e k ) 

kk^E B y = — — r . (id) 

EJVAGN L/i 
i=l 



and 



v 
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As discussed by lEfstathiou et al.l ( 119881 ). a constraint is needed since the likelihood only 
depends on the ratios of the values. We adopt the constraint 



k=l 



- N. 



AGN; 



(14) 



so that the sum of the bins simply equals the number of AGN, and we maximize ln£' = 
ln£ + Xg( £ ), where A is a La grange multiplier . Erro rs in are estimated using the 
information matrix as detailed in lEfstathiou et al.l (119881 ) . In practice, since most of our 
objects have relatively low redd ening values, we pre fer to estimate ^\E BV ) = dn/d\og(E BV + 



0.1). Finally, as discussed by lAssef et al.l (]201lf ). our AGN SED template is as blue as 
possible, so some of the reddening we find is possibly just due to intr insic differences in the 
SEDs of Type 1 AGN. For example, the mean Type 1 SED template of iRichards et al.l (120061 ) 
is similar to our AGN template with E(B — V) ~ 0.05. Since we do not want this to bias 
the results, we subtract 0.05 from all the E(B — V) values before we construct £. 



5.3. Results 



Figure [9] shows the distribution of £'(Ebv) derived using the sample described in §5.11 
The key thing to note is that the distribution falls with increasing reddening, with a dip at 
E(B — V) ~ 2. A smaller, less significant dip may also be present at E(B — V) ~ 0.15. 
Note that there are no objects observed with a best-fit E (B — V) > 14. Luminous objects 



with such high reddening are expected to be rare (see, e.g., lEisenhardt et al.ll2012t IWu et al. 
2012al ; iBridge et al.l 12012k for such extreme cases) . While the general trend of decreasing 



numbers with increasing reddening is in all likeliness real, our sample is small enough that 
the observed dips could in principle be systematic and caused by the non-parametric method 
we used, as it never imposes the requirement of a smooth distribution. However, when the 
sample is divided into three luminosity bins with equal numbers of objects, as shown also 
in Figure [91 the minimum of the dust distribution at E(B — V) ~ 2 appears in all of them, 
further suggesting this feature is real. Assuming that is the case, a few possible explanations 
are possible. 

In the simplest orientation models for AGN unification, most of the obscuration comes 
from the dust torus. However, if the minimum at intermediate E{B — V) is real, it is 
unlikely that the dust forms a continuous medium, as it is very hard to have a physically 
motivated dust distribution that produces such a feature. If the dust is, on the other hand, 
in geometrically and optically thick clouds, the distribution would simply be the distribution 
of the obscuration of the clouds convolved with the distribution of inclination angles and 
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covering fractions. However, this is also unlikely to be consistent with a minimum in the 
distribution at an intermediate obscuration value. Possibly, thick dust clouds are responsible 
for the E[B — V) > 2 obscuration, and these are embedded in a diffuse inter-cloud dust 
medium which is responsible for the lower obscuration part of the £'(Ebv) distribution. 

Alternatively, the two halves of the distribution could be attributed to different sources, 
with the large obscuration coming from thick dust clouds in a torus-like structure surrounding 
the AGN, and the lower obscuration coming from diffuse dust in the host galaxy ISM. Naively, 
one would not expect the distribution of ISM dust obscuration to vary systematically with 
AGN luminosity. When we divide the sample in three bins of luminosity with equal numbers 
of objects, as shown in Figure [91 we observe a significantly different shape for the £'(Ebv) 
distribution in each bin. We consider this as evidence that the dust obscuration is primarily 
coming from the vicinity of the AGN and is hence associated with the torus, and we discuss 
this below in the context of a receding torus. It may be possible that in certain AGN 
feedback scenarios the column density of the residual dust in the ISM left after the AGN has 
gone through the blow-out phase (see, e.g., iHopkins et al.l 120081 ) could be related to AGN 
luminosity during its quasar phase. 

Note that since our sample is inherently magnitude limited, we cannot easily disentangle 
redshift evolution from luminosity evolution. We consider, however, that it is much more 
likely that the evolution in the dust obscuration is primarily driven by the AGN luminosity 
since hardl y any evolution is obs e rved in the UV th rough mid-IR SE Ds of AGN with co smic 



time (e.g., iRichards et al.l 12006c lAssef et al.l 120101 ) . Furthermore, lUeda et al.l (120031 ) has 



shown that the distribution of neutral gas column densities obscuring the X-ray emission of 
AGN is independent of redshift. 

We also investigate the fraction of Type 1 to Type 2 AGN by simply adding up the 
corresponding bins of the £'(Ebv) distribution. We adopt the stan dard X-ray boundary 



of a gas column dens ity of Nh = 10 22 cm -2 (e.g., lUeda et al.l 120031 ) as the dividing line. 
Maiolino et al.l (120011 ) has shown that the value of E(B — V)/Nh is significantly below the 
Galact ic value for AGN, an d varies significantly among different AGN. The median value 
of the iMaiolino et al.l (j200l| ) sample is E(B — V)/N H = 1.5 x 10~ 23 cm 2 mag, which puts 
the Type 1/2 boundary at E(B — V) = 0.15, or A v = 0.47 for R v = 3.1. It also implies 
that our sample does not contain any Compton-thick AGN (Nh > 10 24 cm -2 ). From the 
joint distribution of all AGN, we find that the fraction of objects that would appear as 
Type 1 AGN is 47 ±8%, consistent with an even split between the two types. From a purely 
observational point of view, this is only strictly appropriate for / < 20 AGN given our sample 
selection. However, we note that very little variation in this ratio with /—band magnitude 
is observed in our sample. 
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Figure [TO] shows the fraction of Type 1 AGN as a function of luminosity when we divide 
the sample into three luminosity bins with equal numbers of objects per bin. There is a sharp 
increase in the Type 1 fraction towards higher luminosities. For the lowest luminosity bin 
we find that the fraction of objects appearing as Type 1 is 29 ± 7%, increasing to 46 ± 15% 
for the intermediate luminosity bin, and to 64 ± 13% at the highest luminosity. This trend 
conforms to the idea of a receding torus. Figure [TUI compa r es our observed trends with three 
different models of a receding torus, taken from 
Type 1 AGN by 



Simpson! (120051 ) . We model the fraction of 



fx 



1 



1 + 3 



L 



AGN 



Lag 



N,0 



l-2?/> 



-0.5 



(15) 



which comes from the simple geometry assumed by ISimpsonl (119981 ) and assuming the scale 
height h oc -^agn- At luminosity Lagn,o, AGN are evenly split between Types 1 and 2. 
We first consider the two cases studied by ISimpsonl (120051 ). n amely that of a constant h 
(■0 = 0), and his favored scenario of ip — 0.23. Note that because ISimpsonl (120051 ) used [O ill] 
luminosities as proxies for the accretion disk luminosity, we must fit for Lagn,o, obtaining 
respectively 1.12±g;i| x 10 45 erg s" 1 {ip = 0) and 1.98+°;^ x 10 45 erg s" 1 (tp = 0.23). As 
shown in Figure [10J both of them give a fair representation of the data. If we also fit for the 
dependence of h on Lagn, we find ip = 0.13 ± 0.17 and Lagn,o = l-42jz^°2g x 10 45 erg s _1 . 
Unfortunately, our modest sample size does not allow us to more finely sample the Type 1 
AGN fraction as a function of AGN luminosity and thereby further constrain such models. 
We do note, however, that the reddening distributions shown in Figure [9] have significant 
power to further constrain the dust distribution. This will be be further explored in future 
work. 



Conclusions 



In an earlier study (jPaper I ) we used the extensive spectroscopy and photometry of the 
2 deg 2 COSMOS field to study WISE AGN selection. We found that the simple criterion 
Wl-W2>0.8 and W2<15.05 produces a sample with 95% reliability and recovered 78% of 
the AGN found with Spitzer IRAC imaging to the same flux depth. Here we have extended 
this study using the larger 9 deg 2 NDWFS Bootes field, which has also significantly deeper 
WISE observations than COSMOS. We show that the reliability of a simple color cut quickly 
degrades towards fainter fluxes due to the large number of z ~> 1 galaxies that contaminate 
the color selection. 



Using the extensive UV through mid-IR broad-band photometry available in the ND- 
WFS Bootes field we have studied W2-dependent W1-W2 selection criteria optimized to 
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find AGN at deeper WISE fluxes than those available in the COSMOS field. We provide 
different criteria depending on whether the emphasis is on reliability or completeness. We 
defined a reliability-optimized criteria as a W1-W2 color limit that varies as an exponential 
of W2 2 , where the parameters can be tuned to achieve different reliability levels ( §3.2p . We 
find that for completeness-optimized selection, no dependen ce on W 2 is needed; a simple 



W1-W2 color criterion suffices. We find that the criterion of iPaper II returns samples with 
a completeness of approximately 75%. 

We have also studied the acc uracy of broad-ban d photometric redshifts obtained for the 
WISE AGN candidates using the Assef et al. feoioh basis of low-resolution SED templates 



for AGN and galaxies. We find consistency with the poor accuracy found by previous studies, 
even though our AGN are brighter than those typically used in such studies. Furthermore, 
we find that although the value of the a parameter, the luminosity fraction of the AGN with 
respect to the host plus the AGN, is insensitive to uncertainties in the photometric redshift, 
the best-fit reddening of the AGN component is strongly affected by those uncertainties. 
We have studied the distribution of the best-fit a parameter of the WISE AGN candidates, 
showing they are biased towards high values. This means that WISE AGN selection is biased 
towards objects that are bright with respect to their hosts. S ince the luminosity of t he host is 



roughly correlated with the mass of its central SMBH (e.g.. iMagorrian et al.lll998l ). this can 
probably be expressed as a bias towards AGN radiating at large fractions of their Eddington 
limits. 

Finally, we have studied the distribution of AGN reddening in the WISE AGN can- 
didates. We have shown that although WISE is more sensitive to unobscured objects, it 
still finds considerable numbers of highly obscured objects. Extending the sample to include 
all AGN found over the field with a > 0.5, spectroscopic redshifts from the AGES survey, 
and high S/N WISE W2 fluxes, we have studied the distribution of objects as a function 
of AGN r eddening. We present a formalism based on the step- wise maximum likelihood 



method of lEfstathiou et al.l (119881 ) designed to account for sample incompleteness as a func- 
tion of obscuration. For a subsample of 362 objects with / < 20 and W2<15.73 for which 
the selection function is well understood, we find that the reddening distributions depend 
on AGN bolometric luminosity. The distribution is peaked for unobscured objects and then 
falls relatively monotonically towards E(B — V) ~ 2, raising towards higher values and then 
dropping again towards E(B — V) ~ 10. While it is possible that our small sample size could 
be driving some of the observed structure, we point out that this shape could be explained by 
continuous diffuse dust medium in which optically thick dust clouds are embedded. We find 
that when looking at the complete subsample, 47±8% of AGN are Type 1 (E(B — V) < 0.15; 
see §5.31) . This fraction is a strong function of the AGN bolometric luminosity, consistent 
with the general scenario of a receding torus. At our lowest luminosity bin, centered at 
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-^agn = 3 x 10 44 erg s _1 , we find a Type 1 fraction of 29±7%, which rises to 64±13% for 
the highest luminosity bin centered at Lagn = 4 x 10 45 erg s _1 . Larger samples, such as 
that provided by the combination of SDSS and WISE, will provide greater constraints and 
insight into the dust distribution in AGN. 
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A. Additional Spectroscopic Redshifts in the NDWFS Bootes Field 



The four slitmasks that we observed were designed to target WISE-selected AGN can- 
didates in the Bootes field, though the low source density of such sources allowed for ad- 
ditional spectroscopic targets. We fille d out the m a sks w ith: (1) IRAC-selected AGN can- 
didates, using the two-color criteria of IStern et al.l ( 120051 ) [ Column 1 of table , Targe t Type 
= IRAC AGN]; (2) z > 1 galaxy cluster candidates from lEisenhardt et al.l (120081 ) [Tar- 
get Type = IRAC cluster]; (3) other 4.5 /zm-selected sources from SDWFS, typically se- 
lected to have [3.6] — [4.5] > —0.1 (AB) which efficiently selects galaxies at z > 1.2 (e.g., 



2005 



Kenter et al. 



Galametz et al.ll2012l) [Target Type = IRAC ] ; (4) X-ray sources from XBootes ( Murray et al. 



20051 : Brand et all 120061 ) [Target Type = XBootes]; and (5) MIPS 24 fim 



sources in the field [Target Type = MIPS]. Given the interest and use of the Bootes field by 
a broad community, we include those additional sources here. 

Table [6] presents the results for 129 Bootes sources for which we obtained redshifts, 
including the eleven targeted sources which are also listed in Table The quality flags are 



defined in £14.41 Of particular note is the LRIS mask which confirms cluster 10.220 from the 
catalog of lEisenhardt et al.l (120081 ) to be at z = 0.96. 
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Fig. 1. — Distribution of [5.8]-[8.0] and W1-W2 colo rs in the Bo o tes fie ld. The lines show 
the colors of the galaxy and AGN SED templates of lAssef et al.l (120101 ) . The color of the 
galaxy templates E (red line), Sbc (green line) and Im (magenta line) are shown between 
redshifts (open circle) and 2 (open square), with dots in the tracks in steps of Az = 1. The 
AGN template is shown without reddening (solid blue line) and with E(B — V) = 0.4 (dashed 
blue line), in the redshift interval between z = (open circle) and z = 6 (open square). Dots 
in the AGN color tracks are spaced by Az = 2. The gray dots show all the WISE sources in 
the NDWFS field with W2 < 15.73. 
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Fig. 2.— WISE W1-W2 versus SDWFS [5.8]-[8.0] colors for WISE sources in the NDWFS 
Bootes field. The left panel shows sources with W2 < 15.73, the 10cr WISE detection limit in 
the Bootes field, while the right panel shows sources limited to W2 < 15.05, corresponding to 
the W2 S/N > 10 limit in the COSMOS field. Objects are separated into non-AGN candi- 
dates (light gray dots), WISE and IRAC AGN candidates (blue dots), WISE-only candidates 
(black dots) and IRAC-only candidates (red dots). The median photometric uncertainty for 
each sample is shown in the lower right corner of each panel. 



-35 - 



Completeness/Reliability Fraction 







1 1 1 1 1 J 1 1 1 


0.5 

X 1 1 1 1 1 1 1 1 1 1 1 1 1 


1 1 1 X 1 1 1 1 


1 

1 1 1 1 1 




Y2 (ma, 



Fig. 3. — Reliability (left panel) and completeness (right panel) of AGN candidates defined 
by a > 0.5 selected by a color cut on Wl — W2 as a function of W2 magnitude. Reliability 
and completeness of 90% (75%) are shown as a function of magnitude by the solid (dashed) 
black lines. Objects redder than the top right corner of the panels are missing due to the Wl 
S/N > 3 requirement. The proposed reliability-optimized criteria (eqn.[3]) for 90% (-R90) 
and 75% (-R75) reliability are shown in the left panel by the white solid and dashed lines, 
respectively. The completeness-optimized criteria (eqn.|4]) for 90% (C 90 ) and 75% (C 75 ) 
completeness are shown in the right panel with the same respective line styles. 
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Fig. 4. — Comparison of photometric and spectroscopic redshifts for the i? 90 sample of 
AGN candidates for: a) the full W2 depth sample; b) limited to objects with W2 < 15.73; c) 
further limited to objects with I < 20; and d) even further limited to objects E(B — V) > 0.5. 
Each panel shows the dispersion between the photometric and spectroscopic redshifts for the 
full sample {Az) and for the 95% objects with the best estimates to minimize the effect of 
outliers (Azgs). Black points correspond to objects with |(z p hot — %>ec)/(l + %>cc)| > 0.5. 
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Fig. 5. — Redshift distribution of the i? 90 sample of AGN candidates for the full W2 depth 
of the field (top panel) and limited to objects with W2 < 15.73 (bottom panel). Black 
histograms show objects with spectroscopic redshifts, while gray histograms add objects with 
photometric redshift estimates to the spectroscopic redshift histograms. Shaded histograms 
only include objects with a > 0.5 while open histograms use objects with all a values, 
including the contaminants. 
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Bol / Bol 

Fig. 6. — Distribution of a values for our R 9Q (top) and C90 (bottom) AGN candidate samples. 
Left panels show the full W2 depth samples, while the right panels are limited to objects 
with W2 < 15.73. Black histograms include only objects with spectroscopic redshifts, while 
gray histograms also include objects with photometric redshifts. 
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24680 2468 

AGN Component E(B-V) 

Fig. 7. — Distribution of best-fit AGN reddening E(B — V) values for our R 90 AGN candidate 
samples. Left panel shows the full W2 depth sample, while the right one is limited to objects 
with W2 < 15.73. Black histograms only include objects with spectroscopic redshifts, while 
gray histograms also include objects with photometric redshifts. 
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Fig. 8. — Best-fit SEDs for our sample of highly obscured, high-redshift AGN candidates 
observed with Keck. Only candidates that met one of the selection criteria beyond the very 
inclusive C 90 when using the all-sky data release WISE photometry are shown. 
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Fig. 9. — Distribution of reddening values (solid black circles) obtained after using the non- 
parametric formalism described in §5.2l to account for sample incompleteness. For simplicity, 
we adopt the notation E BV = E(B — V). The top panel shows the distribution obtained 
when using the complete AGN sample, while the lower three panels show the distributions 
in three bins of AGN bolometric luminosity (erg s -1 ) as defined in §2.31 The distribution 



obtained using all objects is repeated as open gray circles in each of the lower three panels 
for comparison. All the distributions are normalized to unity in the lowest Ebv bin. The 
vertical dotted line shows our adopted reddening boundary between Type 1 and Type 2 
AGN. 
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AGN Bolornetric Luminosity (erg/ s) 

Fig. 10. — Type 1 AGN fraction as a function of AGN bolornetric luminosity, as defined in 
§2.31 The solid points are centered at the mean AGN bolornetric luminosity of the bin while 
the open ones are centered at the median value. The luminosity error-bars show the range 
of each luminosity bin. The gray lines shows the best-fit receding torus models described in 
the text. 
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Table 1. Surface Density of AGN Candidates at W2< 17.11 f 



Selection Criterion 


N 


N (a > 0.5) 


Reliability 


Completeness 


Bands 




(deg- 2 ) 


(deg- 2 ) 


Fraction 


Fraction 


Used 


WISE AGN Selection 










R90 


130 


118 


0.90 


0.09 


Wl, W2 


R75 


256 


195 


0.76 


0.17 


Wl, W2 


C90 


3702 


1152 


0.31 


0.89 


Wl, W2 


C75 


2117 


929 


0.44 


0.72 


TTT1 "ITm 

Wl, W2 


Wl-W2>0.8* 


2000 


901 


0.45 


0.70 


ttti Tiro 

Wl, W2 


Jarrett et al. (2011) 


469 


268 


0.57 


0.21 


Wl, W2, WS* 


Mateos et al. (20121 


391 


256 


0.65 


0.20 


Wl, W2, W3 




83 


68 


0.82 


0.05 


Wl, W2, W3, W4 


Asset et al. (20101 


1785 


841 


0.47 


0.65 


Wl, W2 




44 


43 


0.97 


0.03 


Wl, W2, W3, W4 


Wu et al. (2012b) 


3218 


1109 


0.34 


0.86 


Wl, W2 


Other Infrared AGN Selection 








Messias et al. (20121 


662 


543 


0.82 


0.42 


K s , [4.5], [8.0] 




314 


280 


0.89 


0.22 


K s , [4.5], [8.0], MIPS 24/xm 


Stern et al. (20051 


986 


659 


0.67 


0.51 


[3.6], [4.5], [5.8], [8.0] 


Lacv et al. (20041 


2888 


1029 


0.36 


0.79 


[3.6], [4.5], [5.8], [8.0] 


Lacv et al. (20071 


1297 


735 


0.57 


0.57 


[3.6], [4.5], [5.8], [8.0] 



^The effective W2 limits for the Rgo and -R75 criteria are 16.26 and 16.45 mag respectively due to the Wl 
S/N > 3 requirement of our sample. 

*W4 is also used if detected. 



"This criterion corresponds to the color cut proposed by iPaper II without the magnitude limit W2< 15.05 



Note. - WISE AGN selection criteria Rgo, Rj^, Cgg an d C75 are desc r ibed b y eqns. ([3]) and (TJ|. 



The remaining WISE AGN selection criteria are as follows: Ijarrett et al.l (I201lh: W 2 -W3> 2.2, W2 



W3>4.2, Wl-W2>0.1(W2-W3)+0.38, Wl-W2<1.7 and object is not a star Wateos et al.l (I2OI2I) 3-ba nd: 
W2-W3>2.157, Wl-W2>0.315(W2-W3)-0.222, Wl-W2<0.315(W2-W3l+0.796: lMateos et al.l faoi$) 4 - 
band: W3-W4>1.76, W l-W2>0.50(W3-W4) -0.405. Wl-W2<0.50(W3-W4l+0.979: lAssef et al.l (fcoiiih 
2-band: Wl-W2>0 85; Ussef et al.l (|201oh 4-band: W3-W4>2.1, Wl-W2>0.85, W1-W2>1.67(W3- 
W41-3.41; IWu et al.l (|2012bh : Wl-W2>0.57. For the other infrared AGN selection criteria we refer the 
reader to the original studies. 
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Table 2. Surface Density of AGN Candidates at W2<15.73 



Selection Criterion 


N 


N (a > 0.5) 


Reliability 


Completeness 


Bands 


(d 


eg" 2 ) 


(deg- 2 ) 


Fraction 


Fraction 


Used 


WISE AGN Selection 












R90 


117 


106 


0.90 


0.53 


Wl, W2 


R75 


176 


133 


0.76 


0.67 


Wl, W2 


C90 


439 


177 


0.40 


0.88 


Wl, W2 


C75 


194 


144 


0.74 


0.72 


Wl, W2 


Wl-W2>0.8* 


182 


139 


0.77 


0.69 


Wl, W2 


Jarrett et al. (20111 


166 


128 


0.77 


0.64 


Wl, W2, W3* 


Mateos et al. (20121 


161 


129 


0.80 


0.64 


Wl, W2, W3 




69 


56 


0.80 


0.28 


Wl, W2, W3, W4 


Asset et al. (20101 


161 


130 


0.81 


0.65 


Wl, W2 




39 


38 


0.98 


0.19 


Wl, W2, W3, W4 


Wu et al. (2012b1 


347 


170 


0.49 


0.85 


Wl, W2 


Other Infrared AGN Selection 








Messias et al. (20121 


166 


152 


0.91 


0.76 


K s , [4.5], [8.0] 




107 


102 


0.96 


0.51 


K s , [4.5], [8.0], MIPS 24^m 


Stern et al. (20051 


194 


157 


0.81 


0.79 


[3.6], [4.5], [5.8], [8.0] 


Lacv et al. (20041 


484 


183 


0.38 


0.92 


[3.6], [4.5], [5.8], [8.0] 


Lacv et al. (20071 


262 


172 


0.66 


0.86 


[3.6], [4.5], [5.8], [8.0] 



*W4 is also used if detected. 

*This criterion corresponds to the color cut proposed by iPaper 11 without the magnitude limit W2<15.05. 
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Table 3. Surface Density of AGN Candidates at W2< 15.05 



Selection Criterion 


N 


N (a > 0.5) 


Reliability 


Completeness 


Bands 




(deg- 2 ) 


(deg- 2 ) 


Fraction 


Fraction 


Used 


WISE AGN Selection 












GG 


59 


0.90 


0.77 


Wl, W2 


R75 


87 


65 


0.75 


0.84 


Wl, W2 


C90 


120 


69 


0.57 


0.89 


Wl, W2 


C75 


G4 


59 


0.93 


0.77 


Wl, W2 


PaDer I 


G2 


58 


0.94 


0.75 


Wl, W2 


Jarrett et al. (2011) 


GG 


59 


0.90 


0.77 


Wl, W2, W3* 


Mateos et al. (2012) 


65 


GO 


0.92 


0.78 


Wl, W2, W3 




48 


39 


0.82 


0.51 


Wl, W2, W3, W4 


Assef et al. (2010) 


57 


55 


0.96 


0.71 


Wl, W2 




29 


29 


1.00 


0.37 


Wl, W2, W3, W4 


Wu et al. (2012b) 


99 


G7 


0.68 


0.87 


Wl, W2 


Other Infrared AGN Selection 








Messias et al. (2012) 


G7 


63 


0.94 


0.82 


K s , [4.5], [8.0] 




48 


47 


0.97 


0.60 


K s , [4.5], [8.0], MIPS 24fim 


Stern et al. (2005) 


74 


GG 


0.89 


0.85 


[3.6], [4.5], [5.8], [8.0] 


Lacv et al. (2004) 


123 


70 


0.56 


0.90 


[3.6], [4.5], [5.8], [8.0] 


Lacv et al. (2007) 


91 


G8 


0.75 


0.89 


[3.6], [4.5], [5.8], [8.0] 



*W4 is also used if detected. 
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Table 4. Photometric Redshifts 



Sample Az Az 95 bias/(l + z) bias 95% /(l + z ) Nagn z s Fraction 



Full W2 Depth 



R 90 0.27 


0.20 


0.15 


0.14 


618 


0.56 


R 75 0.29 


0.21 


0.14 


0.13 


839 


0.44 


C 90 0.31 


0.23 


0.13 


0.12 


1668 


0.13 


C 75 0.29 


0.23 


0.17 


0.16 


1360 


0.10 


W2<15.73 












R 90 0.27 


0.20 


0.14 


0.13 


595 


0.61 


i? 75 0.29 


0.20 


0.13 


0.12 


731 


0.57 


C 90 0.29 


0.19 


0.09 


0.09 


890 


0.48 


C 75 0.26 


0.20 


0.13 


0.12 


740 


0.48 


W2<15.73 and I<20 










R 90 0.27 


0.20 


0.13 


0.12 


420 


0.85 


R 75 0.30 


0.20 


0.12 


0.11 


507 


0.84 


C 90 0.30 


0.19 


0.09 


0.09 


597 


0.75 


C 75 0.27 


0.20 


0.13 


0.12 


477 


0.82 


W2<15.73 and I<20 and E(B 


- V)>0.5 






R 90 0.25 


0.11 


0.10 


0.11 


68 


0.77 


i? 75 0.42 


0.12 


0.09 


0.10 


90 


0.84 


C 90 0.42 


0.11 


0.04 


0.05 


134 


0.76 


C 75 0.24 


0.12 


0.11 


0.12 


78 


0.78 



Table 5. Summary of Obscured AGN Candidates Observed with Keck/DEIMOS 



IN dlllt. 




R.A. 






Dec. 




I (mag) 


W2 (mag) 


W1 

VV 1 


-W9 
vv z 






d 






E(B- 


-V) 


Z 


n 


/ipf inn 
QClCCllOIl 


INOieS 


W1427+3400 


14: 


:27: 


.54. 


.57 


34 


:00 


:43 


.31 


21. 


.90 


15. 


.81 


0. 


.82 





.386 


± 


0, 


.091 





,41± 


0. 


,11 


1.29 


A 


C75 


[On] 


W1427+3403 


14: 


:27: 


17 


.16 


34 


:03 


:41 


.84 


21. 


.24 


15. 


.76 


0. 


.50 





.209 


± 


0, 


.117 





.33± 


0. 


,12 


1.14 


A 


None 


CaHK 


W1427+3408 


14: 


:27: 


.17. 


.93 


34: 


:08 


:28, 


.60 


21 


.71 


14. 


.99 


2 


.08 





.989 


± 


0, 


.009 





,88± 


0. 


,04 


1.16 


A 


-R90 


[On] 


W1428+3359 


14: 


:28: 


■12. 


,31 


33 


:59 


:25. 


.13 


23 


.22 


15. 


.15 


1 


.02 





.685 


± 





,027 





.87± 


0. 


,05 


1.34 


B 


-R90 


[On] 


W1429+3529 


14: 


:29: 


:54. 


,83 


35 


:29 


:04 


.08 


21. 


.97 


15. 


.77 





.71 





.455 


± 


0. 


,082 





,20± 


0. 


,08 


(1.3) 


F 


C90 




W1430+3525 


14: 


:30: 


:31. 


.69 


35 


:25 


17 


.78 


20. 


.64 


15. 


.07 


1. 


.54 





.944 


± 


0, 


,012 





.45± 


0. 


,06 


1.11 


A 


-R90 


[0 11], CaHK 


W1430+3530 


14: 


:30: 


:00. 


,50 


35 


:30 


:55 


.01 


19 


.42 


15, 


.55 


1 


.06 





.497 


± 


0, 


,044 





,25± 


0, 


,07 


(1.3) 


F 


R75 




W1431+3525 


14: 


:31: 


:06. 


,26 


35 


:25 


:46 


.24 


23 


.63 


15, 


.60 


1. 


.88 


1 


.000 


± 


0, 


,004 





,78± 


0. 


,05 


(1.3) 


C 


R90 


[On] 


W1431+3528 


14: 


:31: 


:31. 


,38 


35 


:28 


:38 


.21 


23 


.62 


15. 


.47 


1, 


,77 





.990 


± 


0, 


,004 





.89± 


0. 


,04 


1.34 


A 


-R90 


[On] 


W1432+3523 


14: 


:32: 


23. 


,02 


35 


:23 


:21 


.41 


19 


.03 


14. 


.33 


1 


.01 





.742 


± 


0, 


,006 





,50± 


0. 


,03 


0.26 


A 


R90 


CaHK,Ha,[Nn] 


W1432+3525 


14: 


:32: 


:37. 


,30 


35 


:25 


12 


.56 


21 


.94 


15. 


,54 


0. 


,72 





.583 


± 


0. 


,036 





.75± 


0. 


,05 


1.12 


A 


C90 


Mgu absn,[On],D4000 


W1432+3526 


14: 


:32: 


:22 


,61 


35 


:26 


:46 


.88 


22 


.85 


15. 


,01 


1, 


,53 





.954 


± 


0, 


,006 





.79± 


0. 


,05 


1.44 


B 


R90 


[On] 



Note. — The AGN selection criteria listed in the last column is the least inclusive one met. For sources that failed to yield spectroscopic redshifts in these 
observations, we list the photometric redshift in parenthesis. Coordinates are J2000. 



Table 6. Additional Results from Keck Observations. 



Target Type R.A. Dec. z Q Slitmask(s) Notes 



IRAC 


14:27:13.39 


+34:09:04.5 


1.343 


A 


C[36] 


IRAC AGN 


14:27:14.32 


+34:09:01.3 


1.692 


B 


C[14l 


MIPS 


14:27:14.63 


+34:08:46.6 


1.343 


B 


C[33] 


MIPS 


14:27:14.77 


+34:06:09.0 


1.151 


A 


C[28] 


MIPS 


14:27:17.40 


+34:07:26.9 


1.082 


A 


C[311 


WISE AGN 


14:27:17.92 


+34:08:28.6 


1.158 


A 


C[03] 


MIPS 


14:27:18.37 


+34:09:01.3 


0.130 


A 


C[34l 


IRAC 


14:27:19.18 


+34:06:47.0 


1.213 


B 


C[37l 


IRAC 


14:27:20.76 


+34:06:42.8 


0.305 


A 


C[38] 


MIPS 


14:27:23.95 


+34:06:15.5 


0.190 


A 


C[29l 


XBootes 


14:27:24.80 


+34:05:12.7 


1.753 


A 


CF151 


IRAC 


14:27:27.55 


+34:05:38.1 


1.483 


A 


C[40l 


MIPS 


14:27:27.56 


+34:04:13.7 


0.126 


A 


C[23] 


IRAC 


14:27:27.90 


+34:06:25.5 


1.690 


B 


CF411 


MIPS 


14:27:29.76 


+34:04:34.8 


1.065 


A 


C[25l 

L "J 


MIPS 


14:27:30.15 


+34:06:05.0 


0.191 


A 


C[27] 


IRAC AGN 


14:27:31.50 


+34:03:28.6 


2.192 


A 


C[13] 


IRAC 


14:27:32.07 


+34:03:19.4 


1.413 


A 


C[42] 


MIPS 


14:27:32.42 


+34:08:08.0 


0.872 


A 


C[32] 


IRAC 


14:27:41.87 


+34:04:55.0 


1.524 


A 


C[43] 


IRAC AGN 


14:27:43.36 


+34:03:06.1 


3.134 


A 


C[12] 


WISE AGN 


14:27:47.17 


+34:03:41.2 


1.137 


A 


C[02] 


serendip 


14:27:47.26 


+34:03:40.8 


1.623 


B 


C[02] 


IRAC AGN 


14:27:47.98 


+34:02:44.1 


3.276 


A 


C[ll] 


IRAC 


14:27:48.28 


+34:03:11.7 


1.427 


A 


C[45] 



Mg II absorption, [O n] 

QSO: Mgn (w/ Mgn absorption system at z=1.342) 
QSO: Mgn 

Mgn absorption, [On],D4000 

[O n], [Ne m] 

[On] 

[O m], Ha 
[O n],D4000 
H/3,[0 in], Ha 
H/3,[0 m], Ha 
QSO: Cm],Mgii 
B2900,[O n] 
H/3,[0 in], Ha 

Mgn absorption, [O II] (do not see past 1/xm in Allslits, to confirm [On]) 

[O n],D4000 

H/3,[0 m], Ha 

QS02: Civ,Heii,Cm] 

[On] 

[Oll],D4000 

Mgn absorption, [On] 

QSO: Lya,Civ,Cin] (self-absorbed) 

CaHK 

[Oil] - v. ft. serendipitous 
QSO: Lya,Civ,Cm] 
[On] 



Table 6 — Continued 



Target Type 


R.A. 


Dec. 


z 


Q 


Slitmask(s) 


Notes 


^A/TCT? a oat 
WlorL AGIN 


1 A .07. K A fin 

14:2 / :o4.dU 


i q /i .nn. /iQ n 
+o4:UU:4o.U 


1.293 


A 

A 


c rm 1 


r<"> ttI 


A TTTDC 

Mil b 


1 A .07. tr cr i o 

14:27:55. 18 


i o a .no. /ii n 

+34:02:41.9 


0.643 


A 

A 


( ~i Tool 

G[22J 


r/^i ttI xj /Q [a tttI 
[U IIJ,Hp,[U IIIJ 


TTD A O 

1KAC 


1 /I .07. r 7 /in 

14:27:57.49 


i o a .no . en o 

+34:1)3:511.8 


n An A 
0.494 


A 

A 


O To rl 
G[35J 


[U IIJ,L-arlK,J4p,ilQ! 


TD A O 

1KAU 


1 A .07. rn £i /I 

14:2 ( :e>9.o4 


i o /i .no .A7 o 

+34:1)3:1) ( .6 


1.621 


A 

A 


\ A nl 
U|49J 


ir\ ttI 


A TTTDC 


14:20:1)3. 55 


+33:59:23.4 


0.756 


A 

A 


ri oi 
C[18J 


TO ttI O , . LJ T. ' T_T Q 

[U IIJ,CaJriK,lip 


MlFb 


14:28:03.65 


i o A no or o 

+34:02:25.3 


n *° Tn 

0.679 


A 

A 


c"\ rr>nl 

G[2UJ 


\ f \ TT "| T T ) \ f \ tttI 

[U Ilj ? Mp , [CJ IIIJ 


IV TTTDO 

MIFb 


1/1 o o n A OA 

14:28:04.84 


i o A no H7 o 

+34:02:27.8 


n /i r* 

0.458 


A 

A 


roi 1 
C[21J 


m ttI tt /3 m tttI tt 
[U Ilj ? Mp , [CJ lllj,lia 




14:28:1)7.7.5 


i 99. ro.oo n 

+33:58:28.9 


1 0£?0 

1.262 


A 

A 


ri 71 
C[17J 


ncn. ai„ tt m TT i tat^ t TT i 
C^bO: Mgll,[U llJ,[JNelllJ 


A TTTDC 


14:28:1)9.59 


+34:01:31.6 


0.627 


A 

A 


Ti nl 

C[19J 


\r\ ttI t_j q ir\ tttI 
[O iiJ,Hp\[<J mj 


A TTTDC 


1 l.oo .nn 71 

14:28:1)9.71 


+33:58:06.6 


1.24b 


A 

A 


C[16J 


rr^ ttI 
[0 nj 


AUootes 


i j.oo.in n7 

14:28:1U.U7 


1 /1 . . 1 n /i 

+34:28:19.4 


1 OO A 

1.234 


A 

A 


b 


a at. / ' tttI r/~* T t1 
AGIN: U IIIJ,[U Ilj 


TTD A O 
IKAL 


i/i.oo.ii no 

14:28:11.28 


1 /I . nn . r /i 

+34:00:54.6 


1 A O A 

1.484 


A 

A 


C'\ rr ol 

C[53J 


\r\ ttI 
[0 nj 


U7TCT7 1 A OAT 


i /I.00.10 on 

14:28:12.31) 


1 . r n . r 

+33:59:25.6 


1 O A O 

1.343 


i ) 

D 


O Tnnl 
G[U(JJ 


\r\ ttI 
[0 nj 


TTD A O 
1KAL 


1 zl.OO.1 /I OC 

14:28:14.2b 


1 A .on. 00 EC 

+34:29:22.5 


0.607 


A 

A 


r 


A OAT 

AGIN 


TTD AO ^,1 , , „j 

1RAO cluster 


i zi.oo.1 /i ni 

14:28:14.91 


1 /i .00.00 
+34:28:28.2 


1 no 

1.08 


D 


r 


t a a nnn 

D4UUU 


IRAC cluster 


1/1 O O 1 C O i 

14:28:15.24 


1 A 07 r 1 

+34:27:35.1 


1 no 

1.023 


B 


T7 1 

1 


[0 11] 


1KAU 


14:28:19.5b 


1 9 /i .09 . 1 n n 
+34:23:10.0 


0.363 


A 

A 


i ■ 

b 


\C\ ttI T_J ^. 

[U HJ,lla 


IRAC cluster 


14:28:21.26 


+34:28:52.3 


0.962 


A 


F 


Mgn absorption, [0 11] ,D4000,H<5 


IRAC cluster 


14:28:22.15 


+34:27:36.4 


1.140 


A 


F 


[0 n],CaHK 


IRAC cluster 


14:28:22.34 


+34:27:21.6 


0.963 


A 


F 


CaHK 


IRAC cluster 


14:28:23.63 


+34:27:17.5 


0.96 


B 


F 


D4000 


WISE AGN 


14:28:25.03 


+34:22:37.2 


1.518 


B 


F 


[On] 


IRAC cluster 


14:28:25.06 


+34:26:59.5 


0.963 


B 


F 


CaK,D4000 


IRAC cluster 


14:28:25.86 


+34:25:33.4 


0.95 


B 


F 


D4000 


XBootes 


14:28:29.06 


+34:25:45.6 


1.036 


A 


F 


AGN: broad Mgll,[Ne v],[0 ll],[0 in] 
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Target Type 


R.A. 


Dec. 


z 


Q 


Slitmask(s) 


Notes 


TD A C 1 


1 A .OQ.Qfl O 1 

14:zo:oU.zl 


I O A .07.90 Q 

+o4:z / :oo.o 


n A A 
U.424 


A 

A 


1 


Mgn absorption, Ha 


IKAU 


1 zi.oo.Qi in 

14:zo:ol.l(J 


i O /I .O /I .OO ^7 

+o4:z4:z5.7 


000 
3.333 


Jj 


1 


asymmetric Lya, Lya break 




1 /I .00.07 OO 

14:z8:o7.oz 


i O /I .O A . C 

+34:24:51). 5 


1 n /i n 

1.049 


A 

A 


1 


[U IIJ,D4U0U 


TD A C 1 


i /i .oo . a n on 

14:zo:4l).o9 


i O/i.O/i .01 n 

+o4:z4:zl.l) 


n ftn /i 

U.694 


A 

A 


1 


AGIN: broad Mg II, narrow lip,|UlllJ 


A/TTT3C 


i a .on. £ £ i a 


+oo:oU:zo.z 


n ftoo 
U.DOZ 


A 

A 


A|4oJ 


[U IIJ,xlp,[U IIIJ (Hp w/ zna component) 


\,r I ) "a. 

Aoootes 


14:29:57.10 


i r on cr /i r 

+35:30:54.5 


0.264 


A 

A 


a ri ci 
A[15J 


/ \ i / \ 1 1 ) r ^ \ , , , 1 tj 

QoU: Hp, [U III], Ha 


serendip 


i /i -on. ftn 
l4.zy.D ( .by 


+00. z 1 ALA 


n /ino 
u.4yz 


A 

A 


A TQOl 

A[oyj 


ro ttI ro tttI 
[U IIJ,[U IIIJ 


Mir b 


14:zy:o 7 .94 


1 9 COT. A O 1 

+o5:z r :4z.l 


n r:/f7 
U.54/ 


A 

A 


a roni 
A|39J 


\r\ ttI tj q rr^ tttI 
[U IIJ,Hp,[U IIIJ 


1RAG 


i /i .on. n /i 

14:29:58.94 


i c .on.07 n 

+35:29:27.9 


1 C7ft 

1.57b 


JJ 


a m 1 
A[51J 


Tf"\ ttI 

[0 iij 


TTJ A C* 

litAG 


i /i .on.nft n7 
14:oU:Ub.U 1 


I QC.Q1 .1 n ft 

+oo:ol:lU.b 


1 o7f; 
l.z/0 


D 


a rcoi 
A[oz| 




Mir b 


14:oU:U f .00 


1 ocon.oQ n 
+oo:zy:zo.U 


1 nn/i 
1.UU4 


A 

A 


A \A Ql 

A[43J 


[U IIJ,JJ4UUU 


TT) A ("1 


14:oU:Uo. <o 


+oo:oU:ol.o 


1 nno 

i.uyz 


A 

A 


a rcoi 


\("\ TT i t ~\ 1 nnn 
[U IIJ,JJ4UUU 


A /TTDC 

Mir b 


i /i .Qn.no on 
14:oU:U9.zU 


I 9 C.OO . K A 9 

+35:28:54.3 


n ncn 
U.U59 


A 

A 


A X A 1 1 

A[41j 


xj q \c\ tttI tj ^, rc tttI 
Mp,[U IIIJ, Ma, |o IIIJ 


TT3 A /~1 A /"* AT 

1RAG AGIN 


i /i .on.nn n £ 

14:30:09.55 


i c .00.00 
+35:28:22.3 


n n^70 

0.978 


A 

A 


A[14J 


ybU-z: lots 01 lines, including [lNevJ342b 


TT3 A C 


i a .Qn-i i no. 
14.oU.ll.Uo 


i or.OV-nQ 1 

+00. z 1 .Uo.l 


1 Cft 


r> 


A[D4J 


Lya (serendipitous?) 


A /TTTJO 

Mlrb 


1 /I .on.i o cn 
14:o0:lz.o0 


+35:2 / :14.b 


n 1 77 
U.l / 7 


A 

A 


A TQftl 

A|3b| 


Hp,[U IIIJ, Ha 


A /TTTDC 


1 /I . on. 1 q 1 /i 
14:oU:lo.l4 


1 9 £ .Oft .CO ^7 

+35:2b:53. 7 


n ofto 
U.362 


A 

A 


A Tool 

A|33J 


uaHK 


MIPS 


14:30:14.33 


+35:26:49.4 


1.481 


A 


A[32] 


AGN: [Nev],[Oii],[Nem] 


MIPS 


14:30:16.54 


+35:30:41.9 


0.792 


A 


A[47] 


[0 n],H/3,[0 in] 


IRAC AGN 


14:30:19.78 


+35:27:04.7 


1.268 


A 


A[13] 


AGN: [Nev]3426,[On] 


MIPS 


14:30:20.75 


+35:26:50.3 


0.681 


A 


A[31] 


CaHK 


MIPS 


14:30:21.18 


+35:29:12.7 


0.0 


A 


A[42] 


M-star 


IRAC 


14:30:25.05 


+35:27:28.8 


0.351 


A 


A [49] 


CaHK,Ha,[Nn] (AGN based on [Nn]/Ha ratio?) 


serendip 


14:30:25.56 


+35:26:28.4 


0.664 


A 


A[29] 


[0 n],H/3,[0 in] 


MIPS 


14:30:26.03 


+35:28:42.4 


0.264 


A 


A[40] 


H/3,[0m],Ha ([On] in 2nd order?) 
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Target Type 


R.A. 


Dec. 


z 


Q 


Slitmask(s) 


Notes 


Mir D 


1 A . on.O^ O A 

14:oU:2o.o4 


i o n .o£? .oo o 
+oo:2b:22.2 


n c\A n 

o.y4o 


A 

A 


A Tonl 
A|29J 


[G IIJjGarlK 


MIFb 


14:3U:2b.9b 


+o5:2b:U5.5 


n o /i n 

0.349 


A 

A 


a ro^7i 
A[27J 


[U llJ,GaHlY,Hp,[U lllj,ria 


NLlr b 


14:oU:o(J. (A 


+35:25:42.8 


1 no /I 
1.U24 


A 

A 


a ro/ii 
A [24] 


rr^ ttI 


Wlbrj AGIN 


1 A .on.oi <y-\ 

14:3U:ol.71 


+35:25:17.8 


i i ne 

1.106 


A 

A 


A Tnnl 

A[00J 


[G llJ,GaHK 


serendip 


i a . on.oi no 

14:o(J:ol.92 


+35:25:17.4 


0.153 


A 

A 


A Tnnl 

A[00J 


ir\ tttI lj ., rc ttI 
[U III], Ha, [b Ilj 


Mir b 


1 A . on. QO QO 

14:oU:o2.o2 


+35:2 / :Ub.2 


1 noe 
1.025 


A 

A 


A|35J 


rr^ ttI n/ir\nn 
[G IIJ,D4000 


A /TTTDC 


14:3U:o7.25 


I O . 0£? . 0<7 1 

+35:26:37.1 


i m o 

1.012 


A 

A 


A Tonl 

A[30J 


/"\o/^\. a /r^. tt ttI xj q \r\ TT t1 
G;bG: broad MgII,[U llJ,J4p,[U Illj 


serendip 


14:oU:o / .zd 


+oo:z / :oo.z 


n eo/i 


A 

A 




\C\ TTl XT ft \f\ TTTl ( ^-^^ ^ v^^-0^ 

[U HJ,rip,[G IIIJ (^ortargetrj 


1RAL 


14:3U:ocj.1U 


+35:24:47.1 


1 one 

1.305 


A 

A 


a red 
A[56J 


\C\ ttI 


Wlbrj AGIN 


i /i .on. A o no 
14:oU:4z.yz 


l Q EC .OQ. 1 A T 

+35:28:14. ( 


1 one 
1.395 


i > 

rs 


A Tnol 
A|02J 


rr^i ttI 


Mir b 


i A - on. /i £ no 
14:oU:4b.Uz 


1 O K .O/l . en /I 

+oo:z4:oU.4 


n ono 


A 

A 


A Tonl 
A|zUJ 


\c\ ttI xj q rr^v tttI xj^ 
[U IIJ,rip,[G IIIJ,rla 


Mir b 


1 A - On. AQ Ad. 

14:oU:4o.4o 


+oo:z / :Uz.4 


n ove 


A 

A 


A To/il 
A [34] 


r^xjxr xj^ 


1KAU AGIN 


i a . on. k 1 nn 


oe.O/l .no n 
+oo:z4:Uo.U 


on/io 


A 

A 


A M 1 1 

A[11J 


t~\Q{^\. r^TlT TTTl A/Trv TT 

t^bG: G IV, U IIIJ,Mg II 


serendip 


i a .on. EC 1 OA 

14:3U:51.o4 


I O CO A .HT o 

+35:24:07.2 


n n 
0.0 


A 

A 


Anil 
A[11J 


star: Ha,CaT absorption 


XT} a a 1\T 

1RAL AGIN 


i /i .on. r 1 ,4 

14:o(J:51.4b 


+35:25:49.4 


n 001 

0.321 


A 

A 


a ri oi 
A[12J 


GaHK 


A /TTTDC 

MIFb 


14:30:56.51 


+35:25:4U.8 


n Ton 

0.720 


A 

A 


A Tool 

A[23J 


[U llj,rsalmer absorption, tip, [U Illj 


A TTDC 

Mir b 


1 /I .01 .nn Qn 
14:31:00.30 


i o e .o e . /i £ n 
+35:25:46.0 


n no /I 
0.0o4 


A 

A 


A roel 
A[25j 


XJ Q \C\ tttI XJ 

rip,[G IIIJ,rla 


MIPS 


14:31:01.91 


+35:24:48.6 


0.650 


A 


A[19] 


AGN: [Nev]3426,[Oii],H/3,[Om] 


MIPS 


14:31:02.47 


+35:24:34.5 


0.651 


A 


A[18] 


[0 n],H/3,[0 m] 


serendip 


14:31:05.86 


+35:24:00.6 


0.414 


A 


A[48] 


H/3,[0 m], Ha 


IRAC 


14:31:06.00 


+35:23:59.4 


0.429 


A 


A[48] 


H/3,[Om] 


WISE AGN 


14:31:31.39 


+35:28:38.3 


1.343 


A 


B[03] 


[On] 


IRAC AGN 


14:31:34.52 


+35:29:23.3 


1.210 


A 


B[15] 


QSO-2: lots of lines 


MIPS 


14:31:34.83 


+35:27:46.6 


0.485 


A 


B[33] 


[On],H/3 


serendip 


14:31:36.26 


+35:29:18.9 


0.163 


A 


B[40] 


H/3,[0 m], Ha 
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Target Type 


R.A. 


Dec. 


z 


Q 


Slitmask(s) 


Notes 


1KAU 


1 A . O 1 . on 7C 

14:31:39. (o 


1 coo.m a 
+35:28:1)1.4 


1.U28 


A 

A 


"D I A Ol 

B[42J 


rr. ttI 

[O IlJ 


Mlrb 


i /i .Qi nn 
14:31:41. UU 


i coo. 1 1 n 
+35:28:11.9 


U.958 


A 

A 


I > fo Cl 

J3|35J 


ro ttI 
[Oil] 


Mil b 


~\ A .Q~\ . A1 1C 

14:31:41.18 


1 con. A A 

+35:29:44.8 


U.163 


A 

A 


H> \ A 1 1 

B[41j 


ir\ tttI rivr ttI rc TT i 
[U lllJ,[JN nj,[b 11J 


sercndip 


1 A .01 - /1 1 OO 

14:31:41. 3z 


1 ocoQ.no a 
+3o:28:U9.b 


1 mi 
1.1U1 


a 

A 


t ) rod 

o|35J 


ttI 
[U IlJ 


Mil b 


14:31:41.8 ( 


l O C •OT'. A n C 

+35:2 r :49.5 


I). 821) 


A 

A 


TO A~\ 

B[34J 


ir\ ttI xj ir\ tttI 
[U IIJ,Hp,[U IIIJ 


Mlrb 


~i A 01 ,4/10 

14:31:44.23 


+35:27:43.3 


0.0 


A 

A 


T> Tool 

B[32J 


star: ual 


MlPb 


1/1 01 Jf7 

14:31:47.43 


1 or 00 n 

+35:28:51.2 


(J. 243 


A 

A 


T> Ton! 

r>[39J 


/ \ , , 1 XT 

[U iiJ,Hq 


A /TTT3C 

Mlrb 


1 A .Ol . /I O 

14:31:48.57 


1 0c.o0.1c c 
+35:28:15.5 


1 01 c 

1.215 


A 

A 


iz? rod 
J3[3b| 


rr~> ttI 
[0 11J 


IRAC 


14:31:52. Ub 


1 c .00.01 
+35:28:21.2 


n c 

(J.7b5 


A 

A 


"D r A Ol 

B[43J 


rr^ ttI r 1 . i_r \c\ tttI 
[U Hj,LaHK,[U IIIJ 


serendip 


1 /1.01 .CO £M 

14:31:53.bl 


1 c o^ .0^; c 
+35:2b:2b.5 


n n 
U.U 


A 

A 


"D I A nl 
r>|49J 


M-star 


lltAL/ AVjlN 


i /i .o i . c c cn 
14:31:55. ol) 


1 c ofi. /i n 
+oo:2b:4U.3 


n 01 1 


A 

A 


"D l"l A~\ 

B[14J 


t^bU: [U HJ,Uarllv,[U IIIJ, broad ria 


serendip 


1 /I.01 .CO rr 

14:31:58.55 


1 0c.oc.c1 a 

+35:25:51.6 


1.136 


A 

A 


B[51J 


ir\ ttI 
[0 11] 


TT3 A f 1 
IRAC 


1 A . 1 - CO O C 

14:31:58.85 


1 c . c . A n n 

+35:25:49.1) 


ncT 

2.957 


A 

A 


13 Id 1 1 

B[51J 


tt. ybU: Lya,L- IV, He II, L IIIJ 


A /TTT3Q 

Mlrb 


1 /I . O 1 .CO 07 

14:31:59.27 


1 O C . /I O 

+35:2b:4b.2 


n o£?n 


A 

A 


I ) foci 

h>[2bj 


[U llJ,Ha; 


TT3 A c 

IKALy 


1 /l .01 . CO C/1 

14:31:59.54 


1 Qc.oQ.01 n 
+oo:zo:zl.y 


n c cv 

U.00 f 


A 

A 


I A A~\ 

r>[44J 


L-axliv 


Mil b 


1 /i . 00 . m 'vo 

14:32:01.73 


i oc.oo J/ 1 n 

+35:28:46.0 


n ono 

l). 392 


A 

A 


T ) fool 

o[38J 


tt \r\ TTT 1 1 1 . 

Hp,[U mJ,Ha 


A /TTTDQ 

Mlrb 


1 /i .00. no no 
14:32:1)3.1)8 


+35:26:46. / 


1). 336 


A 

A 


B[27J 


up ir^t tttI t_i . , 
Hp,[U IIIJ, Ha 


serendip 


14:32:04.97 


+35:26:42.4 


0.438 


A 


B[27] 


[0 n],CaHK,Ha 


MIPS 


14:32:06.77 


+35:25:08.0 


0.732 


A 


B[24] 


[0ii],H/3,[0iii] 


IRAC AGN 


14:32:11.50 


+35:25:34.6 


0.307 


A 


B[12] 


CaHK,H/3,[0 III], Ha 


IRAC 


14:32:13.37 


+35:27:10.5 


0.0 


A 


B[45] 


star: CaT 


IRAC 


14:32:17.71 


+35:25:21.7 


0.368 


A 


B[46] 


[0 n],CaHK,H/3,Ha 


serendip 


14:32:18.14 


+35:25:19.3 


1.136 


A 


B[46] 


[On] 


MIPS 


14:32:19.72 


+35:24:26.5 


0.666 


A 


B[23] 


[0n],H/3,[0 in] 


WISE AGN 


14:32:22.64 


+35:26:46.7 


1.436 


B 


B[02] 


[0 11] (bad slitlet) 
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Target Type 


R.A. 


Dec. 


z 


Q 


Slitmask(s) 


Notes 


WISE AGN 


14:32:23.06 


+35:23:21.4 


0.258 


A 


B[00] 


CaHK,Ha,[Nn] (likely AGN from Ha/[Nn] ratio) 


MIPS 


14:32:26.08 


+35:27:14.7 


0.662 


A 


B[28] 


[On],H/3 


MIPS 


14:32:27.84 


+35:22:42.0 


0.469 


A 


B[20] 


[0 n],H/3,[0 m], Ha 


MIPS 


14:32:28.30 


+35:27:35.1 


0.380 


A 


B[30] 


CaHK,H/3,[0 m],Ha 


MIPS 


14:32:29.30 


+35:22:36.0 


1.011 


A 


B[18] 


[0 n],D4000 


IRAC 


14:32:29.80 


+35:26:29.5 


1.136 


A 


B[48] 


Mgn absorption, [On] 


XBootes 


14:32:32.03 


+35:26:26.6 


1.077 


A 


B[16] 


[0 n],CaHK 


1RAL AGJN 


14:32:36.49 


+35:25:39.4 


1.070 


A 

A 


B[13J 


Q^U: broad Mg II, [U Ilj, broad Banner imes 


WISE AGN 


14:32:37.33 


+35:25:12.6 


1.117 


A 


B[01] 


Mgn absorption, [On], D4000 


MIPS 


14:32:38.38 


+35:21:54.2 


0.495 


A 


B[17] 


[On],H/3,[0 m] 


IRAC AGN 


14:32:39.41 


+35:23:45.0 


1.979 


B 


B[ll] 


QSO: Mgn 


serendip 


14:32:39.81 


+35:21:55.6 


0.835 


B 


B[10] 


[On] 


IRAC AGN 


14:32:40.04 


+35:21:54.5 


0.543 


A 


B[10] 


BL AGN: lots of lines 


MIPS 


14:32:41.29 


+35:26:04.5 


0.071 


A 


B[25] 


H/3,[0 m], Ha 


MIPS 


14:32:45.09 


+35:23:23.4 


0.601 


A 


B[22] 


[On],H/3 



Note. - - Spectroscopic measurements for 129 additional sources. Coordinates shown are in J2000. Q indicates the 
quality of the redshift (see §4.41 for details). Masks A, B and C were observed with DEIMOS, and the bracketed numbers 
indicate the DEIMOS slitlet number. Mask F was observed with LRIS. 



